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XLII. On the Fragmentations of Heavy Cosmic Ray Nuclei 


By K. GorrstsEtn * 


H. H. Wills Physical Laboratory, University of Bristol 
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Max-Planck-Institut fiir Physik, Gottingen + 


[Received December 15, 1953] 


SUMMARY 


In photographic emulsions exposed at 95 000 ft. the reactions of heavy 
nuclei (Z4>3) with the nuclei of the emulsions have been studied. The 
charges and energies of the cosmic ray nuclei were determined both 
‘ directly ’ from their multiple scattering and $-ray densities, or from the 
variation of 5-ray density with range, and ‘ indirectly ’ by measuring the 
opening angles of the showers of fragments, and the total energy released in 
the reactions. A survey is given of the different types of reactions observed. 

The observed frequency of occurrence of heavy fragments (Z>3) 
produced in the encounters of heavier nuclei is in accordance with the 
results obtained by Peters et al. The charge distribution at the top of the 
atmosphere is calculated from the one measured below material (air, glass, 
emulsion) corresponding to 15-40 g/cm? of air, using the collision cross 
section as found by Peters. The influence of energy loss by ionization is 
considered. The resulting charge distribution is in approximate agree- 
ment with the one measured by Dainton, Fowler and Kent, the nuclei Li, 
Be, B being about as frequent as C, N, O, F. The energy distribution 
agrees with the one obtained by Dainton, Fowler and Kent. 

The frequency of meson production in collisions of heavy nuclei is 
studied. 


§1. INTRODUCTION 


THE investigations on the charge spectrum of the primary cosmic radiation 
made by different authors during the past few years by the photographic 
method have led to results not fully in agreement with each other. Thus 
Bradt and Peters (1950) and Peters (1952) came to the conclusion that the 
flux of Li, Be and B nuclei at a geomagnetic latitude of 30° is 10% or less 
of the flux of C and O nuclei, and that the Li, Be and B nuclei observed 
under 15-20 g/cm? of air are mostly fragments of heavier particles having 
suffered collisions inside the atmosphere. Measurements under latitudes 
of 42° and 55° confirmed this result (Kaplon et al. 1952). According to the 
observations of Dainton, Fowler and Kent (1951, 1952), however, the flux 
of the Be and B nuclei at the top of the atmosphere is about 83%, of the 


* Now at the Max-Planck-Institut fiir Physik, Gottingen, Germany. 
+ Communicated by Professor C. F. Powell, F.R.S. 
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flux of GC, N, O and F. The proportion of the Li, Be, B group among the 
primary nuclei is of considerable importance for the conceptions on the 
origin of the cosmic radiation (Bradt and Peters 1950, Dainton et al. 1952, 
Biermann 1953). 

Peters and collaborators determined the charges and energies of particles 
by measuring the grain or 5-ray densities and the residual ranges in stacks 
consisting of alternating layers of nuclear plates and brass or, for high 
energies, the opening and relative scattering angles between the tracks 
of «-particle fragments originating inside the brass layers in collisions of 
their ‘mother nuclei’. Disadvantages of this method are the inaccessi- 
bility to direct observation of events occurring within the brass, and, for 
higher energies, the limitation imposed by the fact that at least 3 «-particles 
must occur as fragments in collisions of primary nuclei for the method to 
be applicable. The Bristol workers obtained their results by measuring 
the grain and 6-ray densities and the mean angles of multiple scattering. 
In consequence they were restricted to tracks of a length of at least 4 mm 
per emulsion. 

The method adopted in the present investigation has already been 
briefly described (Gottstein 1953 a). It consists in the systematic analysis 
of the collisions of heavy nuclei (Z>3) occurring in the emulsion layers of 
nuclear plates. Even if the tracks of heavy nuclei are too short for 
scattering measurements to be feasible, they are very frequently long 
enough to have a number of 5-rays high enough for their density to be 
measured with not too large a statistical error. Considering the numbers, 
masses and energies of the fragments, both of the incoming and the target 
nucleus, and allowing for neutral particles, one obtains the charge and 
energy of the primary in most cases with good accuracy. A further 
advantage of the method is the direct visibility of the events concerned 
thus permitting a detailed study to be made of the phenomena connected 
with the disintegrations. A disadvantage is its restriction to nuclei 
suffering collisions. 


§2. EXPERIMENTAL PROCEDURE 


The emulsions used for this investigation (Ilford G5, 400 uw thick) are 
part of those in which Dainton, Fowler and Kent have carried out their 
measurements. They had been exposed in a vertical position at 95 000 ft. 
at a geomagnetic latitude of 55°. All ‘ stars’ containing ‘ heavy ’ tracks 
with 6-ray densities of at least 1 5-ray per 100 microns have been considered, 
provided these tracks lay in the upper hemisphere and were obviously 
not those of slow secondary particles. Since this 5-ray density corres- 
ponds to the minimum value for Li, all stars being connected with — 
nuclei heavier than He have thus been included. <A ‘ star’ in this con- 
nection was defined as an event in which at least one further track (not 
being a 6-ray) branched off from the heavy track. 118 such events were 
found, the number of tracks per ‘ star’, apart from the heavy primary, 
ranging from 2 to 75. In all there were 2088 tracks in the events. 
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In measuring the 6-ray densities the convention of Dainton et al. was 
adopted, counting all S-rays with 4 or more grains. 

Grain densities, which were measured for all the 1322 ‘ thin’ and ‘ grey ’ 
tracks connected with the events were normalized with regard to the 
minimum grain densities in the respective regions of the emulsions, as 
obtained from the tracks of energetic electron—positron pairs. 

The mean angle of multiple scattering was measured on all tracks of 
heavy primaries and their fragments as well as on all ‘thin’ tracks 
belonging to the events whenever the lengths of the tracks allowed this. 
The coordinate method of Fowler was used (Fowler 1950, Gottstein et al. 
1951) and the influences of spurious scattering as caused by the distri- 
bution of silver grains, mechanical imperfections of the microscope stage, 
distortion of the emulsion and subjective errors were as far as possible 
eliminated by the methods previously described (Menon et al. 1951, 
Gottstein 1953 b). f 

In principle, there exist several methods for the identification of heavy 
nuclei suffering collisions, and as many of them as possible have been 
applied in individual cases. In general, there was good agreement between 
the results rendered by different methods. 

The ‘ direct ’ way consists in the measurement of 6-ray density and 
multiple scattering (Dainton et al. 1951, 1952), or of the variation of 5-ray 
density with the residual range in the stacks of plates (see fig. 4 from 
Dainton et al. 1952). * Indirectly ’, one can obtain the charge and energy 
of the primary nucleus solely by measuring the 5-ray density of its track, 
the total energy of the particle being derived from an analysis of the 
disintegration process. This can be done in one, or several, of the 
following ways : 

(a) Measurement of the Angles 6 between the Tracks of the Fragments and the 
Track of the Primary Particle 


Assuming the mean kinetic energy of protons and a-particles in the 
centre of mass system to be 12 Mev and the greatest energy occurring 
with noticeable intensity to be 30 Mev (Harding et al. 1949) one obtains 


for protons Sato 
V@=—, 
0 
0-25 
ig Onax= aes 
<0 
and for «-particles (Peters 1952) 
ate EOS) 
VA=—, 
0 
0-12 
1 Omnax= ie? ’ 


@... denoting the greatest angle possible in the laboratory system, and «9 
max : : : 
the kinetic energy per nucleon of the heavy primary. To obtain ¢, in 

Bev, 6 is to be measured in radians. 
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Because of the great statistical fluctuations in the case of small numbers 
of fragments this method frequently only gives an upper limit for <p. 


(b) Determination of the Velocities of Fragments 


If there is a narrow ‘ jet’ of particles continuing in the direction of the 
primary, obviously representing the fragments (as e.g., in Plate 8), one 
can in many cases assume that their velocity equals that of the primary. 
If the jet contains particles with grain density between one and four times 
minimum, grain counting alone gives the velocity. In those cases where 
the fragments are in the minimum ionization region without scattering 
measurements being possible, grain counting gives only a lower limit for 
the velocity, but this is sufficient for the identification of the nucleus 
(fig. 4 from Dainton et al. 1952). 


(c) The Smallest Grain Density 


If there is no narrow ‘jet’ of fragments, but there are one or more 
‘thin’ or ‘ grey’ tracks in the event, then the assumption may be used 
that the velocity of the primary was at least as high as that of the fastest 
of the secondary particles (i.e., that with the smallest grain density). In 
general, the fixing of this lower limit is sufficient for the identification of 
. the primary nucleus. The correctness of the above assumption was 
checked in several cases where the primary velocity could be determined 
independently. 


(d) Energy Balance 


The total energy appearing in the reaction must equal the primary 
energy. From fig. 4 given by Dainton et al. (1952) one determines which 
nucleus having the d5-ray density measured, could have supplied the 
required energy. In many cases the result is unambiguous. 

In calculating the total energy of the reaction, a mean energy of 20 Mev 
(including binding energy) was assigned to every ‘ black’ evaporation 
track. The mean energy of a proportion of the ‘ grey ’ tracks was found 
to be 150 Mev (per nucleon) and this value was then used for all ‘ grey ’ 
tracks. ‘Thin’ tracks, whose scattering could not be measured, were 
assigned the energy corresponding to a proton of their grain density. 
‘Minimum’ tracks were assumed to have at least an energy of 1 to 
1-5 Bev. 

To take into account neutral particles, the total energy was assumed to 
equal 2-2 times the sum of the energies of all charged particles. 

In some cases * partial balances ’ were set up, only part of the fragments 
of the primary nucleus being recognizable as such. Apparently one part 
of the fragments had spent their energies for the disintegration of the 
target nucleus, the remaining fragments having lost relatively little 
energy. It was investigated in these cases which assumptions on the 
nature of the individual particles could explain the measurements. 
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§3. RESULTS 
A. General Survey 


As mentioned by Bradt and Peters (1949), three kinds of collisions of 
heavy nuclei may be distinguished phenomenologically : 

(a) The heavy nucleus, having split up the target nucleus, continues 
without loss of charge. 

(6) Several protons and «-particles are split off the nucleus in the 
collision but there is still a heavy remnant (Z>3) continuing approxi- 
mately in the original direction. 

(c) The heavy nucleus is completely split up into protons and «-particles, 
or even merely into protons. 

Plates 7-9 show examples of disintegrations of type (b) and (c).* 

Of the 118 events on which this investigation is based, 2 belonged to 
type (a), 31 to type (6) and 75 to type (c). In 7 cases the primary proved 
to be a relatively slow «-particle, in 2 cases the experimental conditions 
did not allow the determination of the charge of the primary nucleus, and 
in one case the particle having caused the ‘ 8-ray track ’ was found to be 
a fragment of the target nucleus. 

In table 1 are compiled for each primary charge number the fragments 
found, as far as they could be identified. Each compartment corresponds 
to one collision and contains the numbers and types (chemical symbols) of 
fragments identified. The letter ‘H’ in the lower left corner of the 
compartments indicates that there were more than 8 ‘ black’ or * grey’ 
tracks in the events apart from those caused by the fragments, so that 
the target nucleus must have been a ‘ heavy’ one (S, Br, Ag, J). The 
letter ‘ L ’ denotes events with 8 or less ‘ heavy ’ tracks, i.e. collisions with 
H, C, N, O or ‘ glancing ’ collisions with heavier nuclei. The figure in the 
right lower corner gives the energies per nucleon of the primaries in Mev. 
The statistical error in the determination of charges Z>3 was smaller 
than 1 in more than 50° and smaller than 2 in more than 80° of the 
cases. (5 of the 9 nuclei called F could possibly be O or Ne.) In those 
cases where the compartments contain a dash, the fragments of the 
primary could not be distinguished from those of the target nucleus. For 
the lighter nuclei almost all possible modes of disintegration are realized, 
as table 1 shows. There is no significant preponderance of particular 
ways of splitting for any type of nucleus. In view of the variety of target 
nuclei, impact parameters and energies, this is not surprising. 

Of 414 elementary charges carried by the primaries to ‘ H ’ events only 
309 reappear among the fragments, whereas for ‘ L’ events the corres- 
ponding figures are 427 out of 447. (The incoming particles are more 
strongly affected by collisions with heavy target nuclei than by collisions 
with light ones.) 

* The events shown in Plates 7-9 were found in emulsions exposed in a 


balloon flight carried out at Gottingen. Because of the differences in latitude 
and altitude reached they are not included in table 1. 
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Table 2 shows the numbers of different types of identified fragments 
for primary kinetic energies smaller than and equal to or greater than 1 Bev 
per nucleon, separately for H and L events. 


Table 2 
aI : : 
rere Total No. of identified fragments 
H He Li,Be,B CNOF Z>10 
H<1 Bev per nucleon 
H events 27 62 13 9 2 0 
L events 34 39 51 7 5 2 
E>1 Bev per nucleon 
H events 21 100 22 4 0 if 
L events 24 81 36 5 3 0 


As can be seen, the primary nuclei are more thoroughly split up (relatively 
less «-particles and more protons appearing) in collisions with heavy target 
nuclei (H-events) than in those with light ones. This also remains true 
if one considers separately the three groups of primary nuclei with 3<Z<5, 
6<Z<9 and Z>10. (Following the terminology of Bradt and Peters we 
shall call the nuclei belonging to these groups A-, B-, and C-nuclei, respect- 
ively.) The effect is less strongly marked for high energies where even 
collisions with lighter target nuclei result in a more complete break-up 
of the incoming particles. 

One notices further that the relative frequency of occurrence of «- 
particles as fragments is, at least for the L events, dependent on energy. 

Table 3 shows the frequency ratios of «-particles to singly charged 
particles. 

Table 3 
Number of «-particles 
Number of singly charged particles 


A-nuclei B-nuclei C-nuclei 
H events 0-21 0:09 0-31 
L_ events 0-81 1:08 0-49 


These values are in good agreement with those obtained by Perkins 
(1949, 1950 a) in the analysis of the evaporation products of nuclei 
stationary in the emulsion, i.e. ~1 for light nuclei (C, N, O) and 0-4-0-6 
for heavy nuclei, the corresponding ratios in table 3 being 1-08 and 0-49.* 

Figure 1 shows the ratio of elementary charges appearing in ‘ even ’ 
fragments (He, Be, C, .. .) to that in * odd ’ fragments (H, Li, B,...) as a 
function of primary charge for L events and kinetic energies >1 Bev per 
nucleon, these events corresponding roughly, in the opposite frame of 
reference, to the ‘ stars ’ with energetic ‘light’ primaries. Although the 
statistical fluctuations to be expected are rather large, there is a distinct 


ples oe ee 

* Since the collision partners of the nuclei in the reactions investigated by 
Perkins are single nucleons, our L events should be compared with those 
reactions, the light and heavy nuclei of Perkins corresponding to our B and C 


nuclei. 
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tendency for the probability of emission of ‘even’ fragments to be 
higher for ‘even’ nuclei than for ‘odd’ ones. A similar regularity is 
neither recognizable for low energies, nor for H events in general. 


N oven Fig. | 
Nea 
(5) 
oh it 


3° b SEG 7 > 8 ROR 7100 Z 


Number of fragments with even charge number divided by the number of fragments 
with odd charge number as a function of charge number of ‘ mother nucleus’. 
(L-events, kinetic energy per nucleon >1 Bev.) The figures in brackets denote 
the numbers of charged fragments on which the individual experimental points 
are based. 


B. The Charge Distribution of the Heavy Nuclei at the Top of the Atmosphere 


The number of A-, B-, and C-nuclei in a parallel beam having traversed 
h g/cm? of matter is given by the formulae 


Sn he 
W4(4)=W (0) exp (—A/Ag) +N p(0) pa 5 
AA 
(SEBS HAL oer e) tN Cl) a ae 
[exp (—h/A,)—exp (—A/Ac)] 
AES Aso 
+ {3 (0) py pg $28 [ A2e exp (—i/a)—exp (—M/Acl] 
Ap 
— ESE [exp (—h/An)—exp (—h/Ac)] |}, oe ae ae 
ApA 

Np (h)—N a1) ep a hiNe) 7 cea eee = 

x [exp (—h/Ap)—exp (—h/AQ)], Ce SOR rete Oe a ni 


N o(h)=N (0) exp (—hJAc). 1. ey Ue he eee ey 
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N (0), Vz(0), N(0) denote the number of A-, B-, and C-nuclei at the top 
of the atmosphere, A,, Ap, Ag the ‘ absorption mean free paths ’ for the 
different groups of nuclei; jp, is the probability for the occurrence of an 


/ k 
i-type fragment in the break-up of a k-nucleus. 


Ay, ¢ is the ratio of the mean free paths for absorption and collision, for 
B- and C-nuclei 


_ 1 
ee 
uf 


i 


The quantities p, may be derived empirically from table 1. 
k 


It follows, for all events : 


Dx =0-28 Pp, =0-24 
B Total 
Ps =0-26 p » =0-13 
C B+O 
For L events separately : 
Ppa =0-29 p 4, =0-22 
B Total 
Ps =—0-06 Pp » =0-20 
C B+C 


Within the limits of statistical errors these values are in agreement with 
those observed by Bradt and Peters (1950) for collisions in glass (p,—0-23 ; 
B 


pxr=0-24; p 5 =0-23; pono) =0-22) and by Perkins (1950 b) who found 
C B+C B+C 


that in C2 emulsions the average number of heavy fragments (Z>3) 
emitted from stars with 7-11 prongs is 0-20. 

Since the formulae (1), (2) and (3) are valid for a parallel beam, the 
thicknesses of material traversed by each individual particle in air, glass 
and emulsion up to its disintegration have been determined. The mean 
thickness of material traversed by the particles corresponds to 25-7 g/cm? 
of air, if one transforms the individual values of h into those valid for a 
homogeneous medium of air, using the cross sections for the different 
materials as given by Peters (1952). 

In order to obtain the true frequency ratios of the different sorts of 
nuclei, one has to divide their absolute numbers, as found experimentally, 
by their cross sections for collisions. Setting the correction factor for the 
frequency of occurrence of B-nuclei equal to 1, one obtains that the number 
of A-nuclei found should be multiplied by 1-2 and that of the C-nuclei by 
0-75. 

The ‘ projected ’ lengths inside the emulsion layers where the fragment- 
ations occurred have been measured for all primary tracks included in this 
investigation. Table 4 shows the numbers of tracks having different 
lengths. Since the angular distributions of the different types of tracks 
may be assumed to be identical the same should apply to their length 
distributions. It thus follows from table 4 that the relative numbers of 
short tracks which have been overlooked in the scanning of the emulsions 
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are considerably greater for A- than for B- or C-nuclei. We have therefore 
used the assumption that only the values in the last column of table 4 are 
unaffected by this systematic error. 


Table 4 
Projected track lengths (wu) <500 500-3000 >3000 
A-Nuclei 2 13 14 
B-Nuclei 9 25 1h 
O-Nuclei 11 14 6 


Using for A and A the values calculated from the geometrical cross. 
sections and the p; as given by Bradt and Peters (1950) one obtains 
k 


(h=25-7 g/cm?) an abundance ratio of Beryllium and Boron nuclei to 
‘B-nuclei’ equal to 0-79. This is in good agreement with the results of 
Dainton et al. (0-83) whereas according to Peters* the incoming flux of 
A-nuclei (including Li) may be assumed to be smaller than 0-1 times that 
of B-nuclei. 

One obtains further N;,(0)/N .(0)=1-7-+-0-5, a value somewhat smaller 
than that found by the other workers (3-0). 


C. The Energy Distribution of the Heavy Nuclei at the Top of the Atmosphere 


Having determined the thickness of material traversed by each individual 
nucleus before its disintegration, its energy at the top of the atmosphere 
was calculated from the range-energy relation.t Table 5 shows the 
resulting numbers of nuclei in different energy intervals at the top of the 
atmosphere. 


Table 5 


Energy per nucleon 300— 400— 600-— 800— 1000- 


(mev ‘per nucleon) <300 400 600-800 1000 3000 > 3000 
A-Nuclei Oo 4 6 R 2 ff 4 
B-Nuclei O  % 11 53 6 ibe 2, 
C-Nuclei —_—- — 1 5 5 ily 3 


Particles with energies below 300 Mev/nucleon do not occur among the 
A- and B-nuclei; the ‘ geomagnetic cut-off energy’ at the latitude of 
exposure (55°N) is about 350 mev/nucleon (for M=2Z). 

The values of table 5 do not yet represent the true energy distribution at 
the top of the atmosphere because some of the incident nuclei will have 


* As Professor Peters kindly pointed out to us, the value of N4/Ny observed 
by him and his collaborators under 15-20 g/cm? of air was 0-275 and not 0-112, 
. erroneously assumed by us in a previous communication (Gottstein 

953 a). 

} The fact was neglected that some of the nuclei must have been fragments of 
heavier nuclei from collisions inside the atmosphere, but since their number— 
according to §3 B—is not more than 20% of the total, the error involved is 
not serious. 
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lost their entire energy by ionization before reaching the stack of plates. 
To take into account this effect, the distribution of thicknesses of material 
traversed by particles with energies >1 Bev/nucleon has been determined. 
Nuclei with these high energies will not be stopped by ionization loss even 
when having to traverse the greatest thicknesses occurring in our experi- 
mental arrangement. It was then assumed that the distributions of thick- 
nesses traversed by particles of all energies which would have arrived at 
the stack of plates but for ionization loss, were identical. Table 6 shows 
the percentages thus obtained of nuclei of different energies (at the top of 
the atmosphere) which were unable to reach the emulsions because of 
ionization loss. 


Table 6 
Energy per nucleon 300— 400— 600— 800- 1000- ~3000 
(Mev per nucleon) 400 600 800 1000 3000 
A-Nuclei ii 0 0 0 0 0 
B-Nuclei 82 12 0 0 0 0 
C-Nuclei 100 ~—s- 60. il 4 0 0 


If one applies a corresponding correction to table 5 it follows that 9° 
of all nuclei with Z>3 have kinetic energies >3000 Mev/nucleon, and 42% 
kinetic energies >1000 Mev/nucleon. Since ‘ spurious scattering ’ of the 
tracks makes it difficult in some cases to measure accurately scattering 
angles smaller than 0-01 degrees per 100 microns and opening angles of 
fragment showers smaller than 1 degree these proportions may be taken 
to represent lower limits. In table 7 the results are compared with those 
of previous workers, for equal latitudes. 


Table 7 
Proportions of heavy nuclei at the top of the 
Authors atmosphere with kinetic energies 
>1000 Mev/nucleon >3000 Mev/nucleon 
Peters 55% 22% 
Dainton, Fowler, Kent 45% 3%, 
Gottstein 42% 9% 


D. Meson Production in the Collisions of Heavy Nuclei 


The direct determination by multiple scattering and grain-density 
measurements of the masses of singly charged particles is only possible 
if their ionization has not yet reached the ‘ plateau’ value and if their 
tracks within the emulsions are of sufficient lengths. In this manner a 
number of z mesons could be identified among the particles arising from 
the collisions of heavy nuclei. In a few cases particles have been observed 
with mass values similar to those reported by Daniel e¢ al. (1952) for the 
«-mesons, but the magnitude of the statistical errors does not allow 
definite conclusions to be drawn. 

Whenever mass determinations were not possible, the proportion of 
mesons among the ‘ shower ’ particles was estimated in the following way : 
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It was assumed that the number of protons among them is given by the 
difference between the number of charges on the incoming nucleus and the 
total charge of the more than singly-charged fragments, the latter being 
easily identifiable. The remaining ‘ thin ’ tracks were regarded as mesons. 
If the number of the ‘ shower’ particles was equal to, or smaller than the 
number of primary charges, it was assumed that no mesons had been 
produced. 

Of 20 nuclei with energies between 600 and 1000 Mev/nucleon, 6 
produce mesons, 5 of the 6 events being of the ‘L’ type. The average 
number of charged mesons in this energy range is 15/20=0-75. 


Table 8. (Primary energies >1 Bev/nucleon) 


H events L events 
Charge Number of Number of 
Events Mesons Events Mesons 
3, 4 ve 9 4 14 
5, 6 7 19 7 Is 
7,8 3 10 
9, 10 4 19 5 47 
Tak 1b 3 2 Pe 34 
Ip}, 1! 4 4] 1 ll 
15, 16 1 7 1 0 
19, 20 1 4 


For energies >1 Bev/nucleon 33 out of 45 nuclei produce mesons. In 
10 out of the remaining 12 cases the primary charge equals the total 
charge on the singly and multiply charged ‘ relativistic ’ particles. The 
mean number of charged mesons per event for these energies is 


for 21 H events : 4-6 
for 24 L events : 5-5 


There is no marked difference in the figures for H and L events, in agree- 
ment with the predictions of the theory of multiple meson production 
(Heisenberg 1936, 1952). This is borne out by Table 8. 


ACKNOWLEDGMENTS 


The author wishes to express his gratitude to Professor C. F. Powell, 
F.R.S., for the hospitality extended to him in the H. H. Wills Physical 
Laboratory in which this investigation was begun, and to Mr. P. H. 
Fowler for valuable discussions and suggestions. His thanks are also due 
to Dr. B. M. Anand for his assistance with some of the scattering measure- 
ments and 6-ray counts, and to the Bristol team of observers who scanned 
the plates. The microphotographs were made by Mrs. ©. Zollner 
(Gottingen). Drs. M. Ceccarelli and N. Barford have kindly read the 
manuscript. 

It is a pleasure to thank Professors W. Heisenberg and K. Wirtz for 
facilities for continuing this work at the Max-Planck-Institut fiir Physik, 
Gottingen. 


Fragmentations of Heavy Cosmic Ray Nuclei 359 


REFERENCES 


BreRMANN, L., 1953, Annual Review Nuclear Science, Vol. 2. 

Braopt, H. L., and Prrmrs, B., 1949, Phys. Rev., 75, 1779 ; 1950, Ibid., 80, 943. 

Darnton, A. D., Fowisr, P. H., and Kent, D. W., 1951, Phil. Mag., 42, 317; 
1952, Ibid., 43, 729. 

Dante., R., Davies, J. H., Mutvey, J. H., and Perkins, D. H., 1952, Phil. 
Mag., 48, 753. 

Fow er, P. H., 1950, Phil. Mag., 44, 169. 

GoTTsTEIN, K., 1953 a, Naturwiss., 40, 104; 1953 b, Appendix 3 to Kosmische 
Strahlung by W. Heisenberg, 2nd. ed. (Berlin : Springer-Verlag 1953). 

GorttstTEIn, K., Menon, M. G.K., Muuvey, J. H., O’CeaLuaicx, C., and Rocuart, 
O., 1951, Phil. Mag., 42, 708. 

Harpine, J. B., Lattimore, S., and Perxrys, D. H., 1949, Proc. Roy. Soc. 
A, 196, 325. 

HEISENBERG, W., 1936, Zs. f. Phys., 101, 533 ; 1952, Ibid., 133, 65. 

Kapton, M. F., Perers, B., Reynoups, H. L., and Ritson, D. M., 1952, Phys. 
Rev., 85, 295. 

Menon, M. G. K., O’CEatiaicy, C., and Rocuat, O., 1951, Phil. Mag., 42, 932. 

Prrxins, D. H., 1949, Phil. Mag., 40, 601; 1950 a, Lbid., 41, 1388; 1950 b, 
Proc. Roy. Soc. A, 203, 399. 

Prrers, B., 1952, Progress in Cosmic Ray Physics (Amsterdam : North Holland 
Publishing Company). 


[ 360 ] 


XLII. The Ideal Crystal at Absolute Zero 


By L. SatrEer* 
The Clarendon Laboratory, Oxfordt 


[Received December 28, 1953] 


/ 


ABSTRACT 


Theoretical expressions for the internal energy, equilibrium volume, 
and elastic constants of the rare gas solids at absolute zero are derived, 
using a recently proposed approximation for the zero point energy of a 
erystal. The calculations show that, on the basis of a central force 
model, the zero point energy destroys the equality cj.=—c4, (Cauchy 
relation) for these solids. 


§1 
RECENTLY a new expression for the zero point energy of crystals has 
been proposed (Domb and Salter 1952) : 

U,(9/8\NkO, =. 
where @,, the limiting value of the equivalent Debye 90 as T> ~, 
replaces Q,, the limiting value as 7->0. It was there shown that @,, is 
very simply related to the trace of the dynamical matrix for the lattice 
vibration problem, and that for the Bravais cubic lattices in the case 
of nearest and next nearest neighbour interactions eqn. (1) gives quite 
good agreement with the more exact value calculated by the use of the 
moments of the frequency spectrum. 

Assuming, then, that these results may hold more generally (that is, 
for all neighbour interactions and for other lattice structures besides 
the three cubic types), it is of interest to investigate what changes this 
new zero point energy expression implies in the calculations of quantities 
of interest—internal energy, equilibrium volume, elastic constants—at 
the absolute zero of temperature. 

As an example of such calculations we note the work of de Boer and 
Blaisse (1948) with regard to the inert gas solids; they attempted to 
correlate theoretical values with experiment by the addition of a zero 
point energy term U, to the usual lattice energy term Uy). Their 
expression for U, involves the elastic constants (¢,,, Cy), and c,, for the 
cubic crystals) in a rather complicated fashion, so that numerical 
computations are heavy work. There is the further disadvantage that 
the expressions for the elastic constants used are the usual ones (Voigt 
1928, Born 1940) computed from the static lattice energy term U. : only, 
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instead of from the free energy U,+U,. In the following we derive, 
by a unified approach based on the use of the correct free energy, 
expressions for the quantities of interest at absolute zero, including the 
- elastic constants themselves. The resulting equations, simpler in form 
through the use of @,,, confirm and extend the results of de Boer and 
Blaisse. In addition an interesting result is obtained with respect to 
the oft-discussed question of the fulfilment of the Cauchy relations 
(Stakgold 1950). 

We use the formalism of Born (Born 1940) to develop the free energy 
of the crystal at 7’=0 with respect to small homogeneous deformations 
of the lattice and thus obtain, as the coefficients of the strain components, 
the equation of state and the elastic constants. Central forces are 
assumed from the beginning and this allows expression of the results in 
terms of reduced units (de Boer and Blaisse 1948). 


§ 2 
We consider the case of the three Bravais cubic lattices. The lattice 
points of the undeformed lattice are given by r°’=(l,a, J,a, l;a) ; under 
a deformation distorting the cube of side a into a parallelepiped the 
lattice points of the distorted lattice are given by r’=/,a,+/,a,+],as. 
Defining p’ by the equation |r’ |?= |r°’|?+2p’, we can write for p’, 
sure 


tale g 
in Brillouin’s notation. In terms of the a,, e;,=(a; .a,—a6,,)/a”. 


Considering for the moment the unstrained lattice, the equations of 


p'=(a7/2) 1,l,e;;, the e;; being the components of the strain tensor 
1 


eae mi, i= 5S ia, Meare ke cole (2) 
y=1 p=1 3 
where $/,=8,,D¢(|r°"|)+D?4(| r°']), = re a — © bap: 
We can assume solutions of the type 
Die CX (=O) pa oe Oe ee, . (3) 
and also the Born—von Karman plane wave type 
pl=V, exp [i(la)]. exp (—éet), (a)=lyay lary tlyts . (A) 
@,, is given by the expression (Domb and Salter 1952) 
6,,=(h{b[(5/3)ugI?, wo=(BN)Tr[A] 2. (5) 


where 1 is the second moment of the frequency distribution and Tr[A] 
is the trace of the dynamical matrix occurring in the secular equation 
for the frequencies, |4—U?I|=0. The first solution leads directly to 


the result ~ . ; 
Tr[A|= pees ae: SP a ee ate ike \O) 
and the second moment of the frequency distribution is thus 
2a se 
= (120m)? F (x5 +-—)p([r% |. . . . (7 
na=(12ntmyt E (+25) alle) (7) 
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We get, of course, the same result by constructing the 33 matrix [B] 
for the plane wave solution ; in this case p»= (T'r[B])/3, (T7[B]) being 
the value of the trace averaged over all allowable values of the wave 


phases «,. 
The free energy at absolute zero of the unstrained lattice will consist 


of the static lattice energy plus the zero point energy term U,=(9/8)NkO,, 
and by the above relations 


5 1/2 
U,=(3/167)Nk E » | V?d(r) | SS eee oe 
Mi>0 


Sealey 1/2 
and FiN=} J g(r) + — E Be v4 | , =r]... (9) 
I>0 167 Mi>0o0 


We seek an expansion of the free energy per molecule in terms of the 
strain components ; using the expansion (valid for any function f(| r°’|) 
of the lattice sites) 

KV (0? +. 2p) }=f(r) + pDf(r) + (p7/2 D8f(r) + ..--5 +» (10) 


we can easily expand both the quantities ¢(r’) and ( ee 2/r d')r, about 

the undistorted lattice sites r°’. Expanding, in turn, the radical up to 

terms of second order in p’, using the result p’=(a?/2) X1,Le,,, taking 
i) 


into account cubic symmetry, and collecting terms in the e;;, we obtain 


Ff 3h (5S\ V2 
ya} 60+ (=) 


m 


5 \1l2 
jaa {8 (0). = (as) 2B LADglo®) | (exs-beae Fess) 


3h (5S\ 1/2 
(=) [S-1 ¥ 1,4D29(r°") 
1>0 


a 


+ — i6 {Sn + = 
Apa 2’ E 1tDg(r°")}*| (€11?+ €n97+ 33°) 
- 5S\ 1/2 
ete {Sul (=) [S-1 E 1,20,2 D2g(r™) 
m I>0 - 


—1{s- 1 os 1 1° Dg( pol #1} (€44€29-+Coo€as-+ €sae11) 
5 \1/2 
a Ztled'yo®| (€12" + €93°+-31”) +. teres 
| (11) 
Here g(r°)=[3D¢4(r)+1r2D*4(r)],.—,01, S= X g(r), and the S are the static 
I>0 


eres T1804 S ( 


lattice energy terms, as obtained by previous writers (e.g. Born 1940): 
S,M= 3 1,2 Dd(r), 
I>0 


§,,0= a eri 
120 = 21 1212D24(r°). 
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On the other hand, we have the relation from elasticity theory, valid 
for lattice structures of cubic symmetry : 


P/Nv= o— (/2)(€11+€g9+ gs) + (C4/8)(€112 +6 99?-+ C9”) 

+ (Cy9/4)(€11222+ C2033 + 39011) 

+ (C44/2)(€102+-€95?-+ 6332) + 2. hee ee 2) 
Here v=ya® is the volume per molecule, and for the lattices here 
considered the y values are: simple cubic, y=1; face-centred cubic, 
y=2; and body-centred cubic, y=4. f, is the free energy density at 
absolute zero, p is the pressure; ¢,,, C,9, and c,, are the three elastic 
constants for a crystal of cubic symmetry. Dividing eqn. (11) by v and 
comparing coefficients with (12), explicit formulae for the quantities 
of interest at absolute zero are obtained. Adopting the Lennard-Jones 
potential function and writing it in the form 


f(r) = 4e[ (a/7)!2—(a/r)*], 5 icine plies, coe ok 3) 
we have each of these quantities as functions of certain lattice sums 
involving inverse powers of |r°’|. Values of these lattice sums have been 
tabulated (Misra 1940, Barron 1953); we use them in the static lattice 
energy terms, but consider only nearest and next nearest neighbour 
terms in the @,, contributions. Numerically this restriction makes 
little difference (less than 4°/), since terms arising from the dynamical 
matrix are insensitive to more than the first few rings of neighbours. 
We use reduced units (the reducing parameters « and o appearing in the 
force law (13)) to express the results, writing for the case of the face- 
centred cubic structure V*=(v/o), T*=(kT/e), p*=(o%p/€), A*=h/ov/(me), 
C,;*= (o%c,,/€), a=o(V*/2)1/8, and obtain, finally, the following expressions : 


fo 6066 14-454 2-291A* 


Ot = a = pen — pee pars mea pene =e ie 2 + ae (14) 
4 DR. 2299 A* 5 a 5 *2 

meron | 038221 (ee er fee E15) 
Ve V3 V*10/3 Riz 

ae 171-2 123-0 | 1-146A* 1,7. 
V5 Vx3 V*10/3 

19-40—3-320V *? 1 /59-96—16-0V*? E (16) 
x | 8275 SS srs R ee se es 
84:25 54-11 1:1461* 51/2 

ieee Ves, a pease ~ elisa: 

9-625— 1-563 V *? 1 /59:96—16-0V*? ll (17) 
[saa ae PR Rothe 
84:95 54-11  9:483A*(9-625—1-563V*?) _ (18) 

Cee es Fe V¥10/3 Ri ; 
0:06367A* (59-96—16-0V*?)? _ 
Ca*—Cyo ~~ Ta 510s abe ee es (19) 
where R=4-283—2-00V *2, 


Z2B2 
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The use of the reduced units (de Boer 1940), where the parameters € and o 
are experimentally determinable from gaseous phase measurements, 
facilitate comparison with experimental results and in particular with 
the previous theoretical work of de Boer and Blaisse (1948). 


Fig. 1 


§3 
Figures 1 and 2 are after de Boer and Blaisse. Equation (15), with 
p*=0, becomes a relation between V*, the equilibrium volume at Zero 
temperature, and the parameter A*. Corresponding values are shown 
in fig. 1 as the smooth curve, while the points represent experimental 
values for the inert gas solids. A* may be interpreted as the reduced 
de Broglie wavelength of relative motion for a substance fulfilling certain 
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restrictions, and is a measure of the influence of specifically quantum 
effects on the behaviour of the substance. When corresponding values 
of V)* and A* so obtained are substituted in eqn. (14), the curve for U)* 
the latent heat of vaporization at absolute zero, as a function of A* fg 
the result as shown in fig. 2. 

There is reasonable agreement to within the neighbourhood of 
the A* value for deuterium, but the calculations break down for 
larger A* values. The Vy*-A* curve passes through a maximum A* 
value of 1-22 and then curves back rapidly toward the V>* axis, giving 
A*=0 at V>*=1-46. This behaviour is due to the variation of the factor 
(4:283—2-0 V**)1/2 (essentially the volume dependence of @,,*); for 
reduced volumes greater than 1-46 @,,* is imaginary. This is the same 
sort of difficulty which occurred in the calculations of de Boer and 
Blaisse, where a @, computed from theoretical values of elastic constants 
was used for the zero point energy. One would expect, in the light 
of recent developments in the theory of the ideal crystal (Born 1951) 
that such difficulties are inherent in any approach based on the 
conventional ‘ semi-classical ’ lattice dynamics ; the zero point vibrations 
are so large for the solid elements with a very high A* value (H,, 4He, *He) 
that the dynamical assumption of small vibrations may be no longer valid. 

The comparison of theoretical with experimental values of V)* and U)* 
is not a very sensitive test for a zero point energy expression, partly due 
to the uncertainty in the experimental values themselves. The agreement 
obtained using 9,, seems to be as good as that obtained by any previous 
calculations, and has the advantage of simplicity insofar as the actual 
numerical calculations are concerned. 


§4 

Of greater interest, perhaps, is the fact that the present theory yields 
an explicit expression for the elastic constants at absolute zero, consisting 
of the usual ‘ static’ term, as obtained by Voigt and many others, plus 
an additional contribution due to the zero point energy term occurring 
in the free energy. This additional contribution is such, assuming that 
the zero point energy is correctly taken into account by the use of @,,, 
as to destroy the validity of the Cauchy relation cj,=c,4 for the cubic 
lattices. Theoretical calculations of the equation of state and elastic 
constants at high temperatures have been made for simple monatomic 
substances crystallizing in the body-centred cubic (Born 1939) and 
face-centred cubic (Bradburn 1943) structures ; it was emphasized that 
for any finite temperature 7' the Cauchy relations no longer necessarily 
hold. It may not, however, be generally recognized that at absolute 
zero the presence of zero point energy may be sufficient to vitiate the 
relations. 

In fig. 3 theoretical values of the reduced elastic constant Cy, for the 
face-centred cubic lattice are plotted as functions of reduced pressure p* 
for various values of A*. The A* values correspond to the inert gas 
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solids A, Kr, Ne and D (dotted curves). The relative difference 
(C44 —C49*)/C19*=(C44—C19)/C1g is plotted against pressure on the same 
figure for the same elements (solid curves). For example, the theory 
predicts that for the element Kr, at 0°K (physically, at temperatures 
low enough that the magnitude of contributions from the (7/@)* and 
succeeding terms in the free energy expression are negligible), c,, has 
a value of (3-31)(101°) dyne/em?, while ¢44=(3-54)(101°) dyne/cm?, a 
relative deviation of approximately 7%. For Ne, at zero temperature 


Fig. 3 


200f 


and a pressure of 450 atmospheres, the relative deviation is about 
54%, with c,,.=(0-916)(101°) dyne/em?. Provided the elastic constant 
measurements are feasible, such differences should be detectable. The 
experimental difficulties are great, since the usual elasticity techniques 
must be combined with low temperature technique ; there seem to be 
at present no measurements available for the substances to which the 
present theory is applicable (but see Barker et al. 1953). It is of interest 
to point out that in experiments on solid helium carried out at the 
Clarendon Laboratory Dugdale and Simon (1953) sueceeded in determining 
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the bulk modulus K=(c,,+2c,,)/3 as a function of temperature and 
volume, and found a constant value at constant volume for temperatures 
such that 7/0<0-09. In other words, for solid helium the range of 
temperatures such that only the static lattice energy and zero point 
energy terms of the free energy are appreciable is accessible experimentally. 
The solid elements Kr, A, Ne, etc. have @ values greater than that of He ; 
it is thus conceivable that the predictions of this theory are capable of 
verification by measurements carried out at moderately low temperatures. 


§5 
Since the zero point energy makes a contribution to the magnitude 
of the elastic constants a deviation of the measured value from the usual 
theoretical value (calculated from the static lattice energy only) affords 
a quantitative test for any zero point energy expression (Simon 1934, 
Megaw and Simon 1936, Megaw 1939). Considering in particular the 
bulk modulus, for it is the most easily measured, we write 


K=K,+K,, . . . . . . . . (20) 


K,, being the static lattice energy part and K, the zero point energy part. 
We note that at absolute zero the isothermal and adiabatic moduli of 
elasticity coincide, so that there is no need to distinguish between them 
here. For the face-centred cubic lattice, in reduced units : 


at AG RIGO ABV 


kee Sama Ries Te (21) 
1-146A* 106-6— 17:83 V*? 1 /59-96—16-0V*?\2 
i V*10/3 Riel R - a ( R 3 
eetOnses 2-0 VF ke Slane yr irr wh (99) 


For argon at a pressure of 1800 atmospheres: K,/(K,+K,)=0-076 ; 
for neon at 5650 atmospheres: K,/(K,+K,)=0-52. No numerical 
values are given for K, and K,, since such values are sensitive to the 
form of the theoretical potential law adopted. For example, the potential 
o(7)—=ae—" —br-*is preferable on theoretical grounds to the Lennard-Jones 
force law and may be expected to lead to somewhat different magnitudes 
for the c,; and K, though reproducing quantitatively the proportional 
influence of static lattice and zero point energy terms. Pending a further 
investigation based on the use of different potential laws, we have here 
limited the quantitative considerations to the relative effect of the two 
energy contributions. 
Provided, then, that experiments on the bulk modulus can be carried 
out in the ‘ plateau’ region, where K is constant before the temperature 
dependent terms become effective and the value begins to drop with 
rising temperature, measurements on the lighter inert gases at medium 
pressures should give important information on the question of Zero 
point energy and its influences on the elastic properties of the solid. 
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The question of the validity of the Cauchy relations is important in 
the following connection: one of the usual criteria for the existence of 
central forces is the fulfilment of the Cauchy relations. The experimental 
fact that in metals the relations are generally not fulfilled is taken as 
confirmation of the supposition that metallic binding forces are chiefly 
non-central, a supposition which is of course amply substantiated by 
theoretical considerations. Any deviation due to the presence of zero 
point energy one would expect to be negligible in the case of metals, 
for experiment indicates that with the exception of the lightest solids 
the influence of zero point energy on the properties is negligible. It is, 
however, precisely in the case of these exceptional, inert gas type solids 
that the assumption of central forces has the strongest theoretical basis. 
Should experiment show that for these ‘ideal’ monatomic: solids the 
Cauchy relations are not valid, this would not necessarily invalidate the 
central force hypothesis; the present calculations are based on such 
an assumption and yet indicate that in general such a deviation could 
be expected due to the influence of zero point energy. 
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ABSTRACT 

The Simon melting equation p,,/a=(T,,,/T'y)°—1 is derived assuming the 
validity of the Griineisen equation of state and the Lindemann melting 
formula along the melting line. A relation between c and Griineisen’s + 
is derived and discussed in view of available experimental data. Poor 
agreement is found in the case of the alkali metals, the agreement for 
argon being satisfactory. A consequence of possible geophysical interest 
arising from the theory is discussed. The detailed data available for 
solid helium are used to study the assumptions underlying the theory. 


$1 
THE purpose of this paper is to point out a relation between three well 
known equations connected with the solid state : 
The Griineisen equation of state 


aus y 
p= AVanvs ae 0 (1) 
the Lindemann melting formula 
mV 2/362 
Ue eae we eer ee eee” ots) (2) 
and the Simon melting equation 
Pn (em \ 3 
2 = (7) Peat aR Tit =. fede! (3) 


In (1) p is the pressure, U, is the internal energy at absolute zero, 
U—U | C,,dT is the thermal energy, 6 (not necessarily the Debye or 
0 


Einstein 0) is a function 6(V) of volume only such that C,=f(T/@), and 
y=—din6/dinV is Griineisen’s ‘constant’. In (2), @ is again a 
characteristic temperature of some kind, V,, and 7',, are simultaneous 
volume and temperature values along the melting line, and m is the 
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atomic mass. In the Simon equation (3) p,, and 7',, are the melting 
curve values of pressure and temperature, 7’, is the triple point tem- 
perature, and a (essentially the internal pressure) and ¢ are constants. 

Now egn. (1), though originally derived by Griineisen in an extension 
of the Mie theory of solids and based on a particular model, is in a sense 
a thermodynamic relation. That is, given the relation C,=f/(7/6), 
with @ a function of volume only, eqn. (1) and the equally important 
relation 


eee 4 
a 8 ee eed 


(« the compressibility at constant temperature, « the volume expansion 
coefficient) follow at once by simple thermodynamic reasoning (see, for 
example, Davies 1952). Theoretically, the restrictive condition on C, for 
a Griineisen solid implies that the 3N eigenfrequencies (for a monatomic 
solid) of the dynamical matrix are independent of temperature and vary 
identically with volume. Griineisen (1926) has shown that experiments on 
a number of solids are adequately described by the theory. (It has recently 
been shown by Bijl and Pullan (in print Phil. Mag.) that Grineisen’s 
relation breaks down at temperatures below @; we are here, however, 
concerned mainly with much higher temperatures only.) 

The Lindemann equation (2), while originally derived (Lindemann 
1910) on the basis of a simple model of the solid and a particular picture of 
the melting process, is undoubtedly of wider validity ; it finds justification 
in that a number of substances are found to satisfy the formula at the 
normal melting point, while recently Dugdale and Simon (1953) have 
found that solid helium satisfied eqn. (2) over a wide portion of the 
melting curve. Values of the constant vary somewhat from one type of 
solid to another, but for the majority of solids the constants lie within a 
small range of variation. 

The melting curve eqn. (3) was proposed by Simon (1929) after a 
careful study of experimental data; it is found that for all substances 
heretofore studied an a and a c can be determined which will make 
In (p,,+-@) and In 7’, vary linearly. The equation has been very useful as 
a means of representing experimental results for a wide variety of sub- 
stances over varying melting temperature and pressure ranges, and in the 
planning of further experiments (e.g., the use of solid helium as a ‘ model ’ 
substance (Simon 1937, 1952 a)). It is somewhat surprising, therefore, 
that a relation of such wide experimental validity should have, apparently, 
little theoretical foundation. Of the few theoretical derivations which have 
been published, that of Domb (1951) as improved by de Boer (1952) is 
the most successful ; it is based on an adaptation of the Lennard-Jones 
and Devonshire (1939) model of the melting process. The results are 
restricted in validity by the assumption of high pressures and ‘ classical ’ 
substances (negligible zero point energy, 7',,/0> 1). 

We shall show that eqn. (3) follows in an elementary fashion from (1) 
and (2) on the basis of certain reasonable assumptions. 
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First, we assume that the Griineisen equation of state holds along the 
melting curve. We write, then, for (Lye E 


Pm= =U ( Vin) ? lies Uo) rns Vin’ sf) 


Obviously, by eliminating V,, between (5) and the Lindemann equation 
(2) we will be led to a melting curve equation of some form. It is first 
necessary to consider the physical significance of the parameter a occurring 
in the Simon relation. Simon (1952 b) has put forth strong evidence that 
a is related to the internal pressure (9U/AV), of the solid ; in particular, 
a=(0U/0V)p_». For, writing eqn. (3) at 7,0 (assuming that this 


extrapolation is valid) pO. Se ee pe (6) 
On the other hand, we have the thermodynamic relation 
0U 
(57),=2 (Sn), -” Poet o eee ooe (7) 
and in particular, on the melting curve intercept at 7',,—0 
0U 
(57), -o =P et ws ee tom a8 (8) 
Comparing (6) and (8) we are led to the relation 
yay ples 
Chee Deane 


We apply these considerations to the Griineisen equation (5) ; by definition 
(QU/0V)p-p=U, (V), thus 

PEED. hee vo oles a) 
We assume in all that follows that U,'(V,,) is constant along that portion 
of the melting line corresponding to low melting pressures. If melting in 
the vicinity of the triple point (essentially atmospheric pressure) takes 
place at volumes near the inflection point on the U,(V) curve, this should 
be a reasonable assumption. Naturally, at high melting pressures the 
volume will be decreased and U,'(V,,) will vary, but this variation in the 
theoretical a would be completely masked by the large value of p,,. We 
note that the identification (11) was the essential improvement made by 
de Boer (1952) on Domb’s derivation of (3). Now eqn. (5) becomes 


Pm ta= g-(U-y). 2 2  () 


m 
To proceed further we must assume some definite functional dependence 
of U—U, on J, and T,,, ; the limiting high temperature expression 3Nk7',, 
should suffice for most solids (we have in mind throughout monatomic 
solids, for the extension is obvious), the exceptional case of solid helium 
being reserved for later discussion. If we now assume that y is constant, 


Bo that’ we may write ~ 9 p(hikmiA\V Yen... ad « 0 ee woth) 
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b a constant determined by the lattice structure and the interatomic 
potential, then, substituting in (11) and eliminating V,,, between 


PntO=3NETy|Vpqs 8=0(h|kyt!®)V~, 
and m V,,,,2/3 62/T',,=h?.B/k (B dimensionless), we find 


Bk\e-1 
pn ta=3Nky (Fe) ToS nL 3) 
6y+1 
=s— > mnier an acral les 
with ty $2: (14) 


Applying this at the triple point 7'9, where p,, may be neglected in com- 


parison with a, Bivet 

a=3Nky (=) fh mt tone es, 
and dividing (13) by (15) we establish p,,/a=(T,,,/T'))°—1, the Simon 
formula. 

This derivation, then, connects the melting curve parameter c with 
Griineisen’s y, the latter being determined by the form of the interatomic 
potential energy function and the dynamical behaviour of the lattice. We 
note that it has been assumed implicitly in the foregoing, in addition to the 
constancy of y and U,(V,,,) along the melting curve, that the character- 
istic temperature 6 occurring in (1) and (2) is the same function, at least 
in so far as volume dependence is concerned. 

By means of the Clausius—Clapeyron equation we obtain, differentiating 
(3) and using (15) 


Dm we { Aa \ eel a ore — (6by+1\ (BkT,,\e-1 
a= (at) = 7, (re) =28y(G5) Fat) «a9 
and at the triple point 7', 
dPmn\) (Ada \ > ass, (OV ted \ a) Be aN a 
(az), = avs), 28 (Bs) (Feet). an 


Equations (16) and (17) give further relations involving the y of the solid 
phase and experimentally determinate quantities associated with melting. 


$3 

Professor Simon has pointed out to the writer that the relation (14) 
connecting c and y is of interest in the following regard. The condition 
that the entropy should increase along the melting curve with rising 
pressure, or alternately, that melting of a substance by the application 
of adiabatic compression alone be impossible, may be reduced to the 
demand that 7’,,/0 should increase with p,, (Dugdale and Simon 1953). 
The present theory enables us to express this condition in terms of y 
for the substance. For 


a Py 2? 'n 7" ed f if 
dp, \ 2) tae ar (#) 


= 3by—2 1 T' 4 fA 
Sh y=) eas (=) 20 for y222/3. se on 18m 
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The limiting value y= 2/3 is, somewhat artificially, the limiting theoretical 
value for V — oo (ie., the perfect gas), as has been pointed out to the 
writer by Dr. C. Domb. Here y for a perfect gas is defined by the 
relation (4), or equivalently, by y=—(@ In 7/@ In V)s. Using (14), the 
restriction y>2/3 may be expressed as c<5/2. The only metals for which 
experimental c and y values are available are the alkali metals ; without 
exception the experimental values available satisfy y>2/3 but c>5/2. 
As stated below, there is some uncertainty in the y values, so that the 
implied inconsistency in eqns. (14) and (18) may be only apparent. 
Pending further investigation of the question, it may be pointed out 
that the entropy of solid sodium has been shown from experimental data 
to increase along the melting curve (Simon 1937). The case of greatest 
interest geophysically is iron ; Griineisen (1926) found by calculation, using 
eqn. (4), y=1-60, so that the condition for entropy increase along the 
melting line is satisfied. 


$4 

Broadly speaking, a relation between y and c of the form (14) does seem 
plausible in view of available experimental data. Metals, which have 
low y values, have values of ¢ ranging from 2 to 4 (but again, we have 
melting curves for only the most compressible metals). The inert gas 
solids, on the other hand, have relatively high y’s (if we take the lone 
value for argon as representative) and c=1:5. Quantitative checking of 
the theory is made difficult by the fact that for most substances y must be 
computed indirectly from the experimental data; this may be done in 
several different ways and the values so obtained are not always in good 
agreement. At the same time, experimental difficulties in some cases 
preclude melting curve measurements of sufficient accuracy to determine 
unique values for c and a. 

In table 1 are given experimental values of y, the theoretical values 
of c computed using eqn. (14), and the experimental values of c. 

The agreement in the case of the alkali metals is distinctly poor, and 
we conclude that for these highly compressible metals the assumptions 
underlying the present theory are not fulfilled. N evertheless there seems 
to be a correlation in relative magnitude between experimental and 
theoretical c’s (only Na is out of line). We have used Griineisen’s 
calculated values of y, computed through eqn. (4), as being the most 
reliable in the absence of direct determinations. ‘These depend critically 
on compressibility values ; in fact, it may well be that for these metals 
the theoretical calculations of y must be altered in accord with finite 

i eories (Truesdell 1952). 
Saar fe Tea that the substances most likely to fulfill the 
assumptions would be the heavier inert gas solids : A, Kr, and Xe. 
Unfortunately, the only y value available (with the exception of the direct 
measurements for solid He as discussed below) is for argon, as listed in 
table 1, The good agreement between computed and theoretical values 
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of c in this case is blighted somewhat in that the value y=4-5 is suspect 
as being too high (private communication from Professor Simon). 
However, correction to the more probable value 3-5 only betters the 
agreement, giving c=1-16 in agreement with Robinson’s measured value. 
Experimental y measurements for these rare gas solids would be of 
great interest, not only in checking the present theory, but because y 
as a function of V is very sensitive to the form of the interatomic potential 


energy curve. 


Table 1 

Substance y Ctheor Cexptl 
Na 1-25 (1) 1-55 3:56 (7) 
K 1-34 (1) 1-50 4:53 (7) 
Rb 1-48 (1) 1-44 4:2 (7) 
Cs 1-29 (1) 1-52 4:75 (2) 
A 4:5 (3) 1-12 1:16 (4) 

1-228 (5) 

1-478 (6) 

References : Griineisen, 1926, Handb. der Physik, 10. 


(1 

a Simon and Glatzel, 1929, Zeits. anorg. allg. Chem., Bd. 178, 309. 

(3) Simon and Kippert, 1928, Z. f. phys. Chemie, 135, 113. 

(4) Robinson, 1952, Thesis, Oxford. 

(5) Simon, Ruhemann and Edwards, 1930, Z. f. phys. Chemie, B, 
6, 331. 

(6) Bridgman, 1935, Proc. Am. Acad. Arts and Sciences, 70, 1. 

Note: Of the three c values listed for argon, that of Robinson is 

probably the most accurate, as his experiments covered a wider 

melting pressure range. 

(7) Data by Bridgman, 1949, The Physics of High Pressures, 

London. c values computed by Dr. J. A. W. Huggill. The 

values for lithium have been omitted as only a rather 

impure sample was available. 
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§5 

Equation (14) gives a relation between a and 7’, for a particular solid, 
involving the constants c, b and B. On theoretical grounds (de Boer 
1948) we can expect that 6 and B depend only on the lattice structure 
and on the shape of the potential energy curve; by ‘shape’ is meant 
that the potential function can be written in the form 4(r)=eg(r/c), where 
g(x) is a universal function specifying the form of the potential, « and 
a being reducing parameters characteristic of the particular solid which 
specify the scale. The properties of the heavier inert gas solids seem to 
be well described theoretically through the use of such a ‘reduced’ 
potential function ; it is known that they crystallize in the same structure, 
the face centered cubic. Re-examining our equations in the light of 
these dimensional considerations we write for the constant b occurring 
in eqn. (12) . Be 


1/2 
(No®)”, 


oO 
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where now £ is dimensionless and should, for corresponding solids, 
depend only upon the lattice structure. We can express the equations 
in reduced form, using the variables V*=V/No*, p*=po%le, T*=kTle, 
and A*=h/oy/(me) (de Boer 1948) and derive as before 


B e-1 
Oy 4 0.37 (ay) (eeu tas . . . . (13) 
with c=(6y+1)/(6y—2) as before. Equation (15) becomes 
B e-1 
a*=3y (saya) Wea Gr. thecy bags (15’) 


Now Clusius (1936) has determined the constant B for A, Kr, and Xe, 
finding the almost identical values (1-39)(1017), (1-40)(1017), and (1-41)(1017) 
respectively. Since these solids have the same lattice structure, we 
expect the f’s and y’s to be identical. Hence we can write (15’) in the 


f 
Bales PESTACIESC «TO EY OE aay) 


K and c the same constants for each of the corresponding solids. These 
theoretical considerations lead us to expect that experimental values 
of c will be very nearly the same for these three rare gas solids, and 
(using Robinson’s a and ¢ values for argon, known triple point tempera- 
tures, and values of « and o determined from gaseous phase experiments) 
that a for Kr and Xe will be approximately 3900 and 4400 atmospheres 
respectively. Clusius (1940) has determined the melting curves for Kr 
and Xe, but over too small a pressure range to allow a unique deter- 
mination of a. 

Exactly the same considerations apply to eqn. (17); in reduced form 


AS* 1 (AH* Bee! 1 
(ar+) ,.— 7 (as) = (am) (= 17 


1 AH* ja OO 
CRAG Ci eee (a) ==CONSt. es. 4 (20) 


Experimental values of 4H and AV are available (Guggenheim 1950) 
and, using Robinson’s value of c for argon, we get the following computed 
values of (7',*)-*(4H*/AV*)p» for A, Kr, and Xe: 12:0, 11-7, and 
11-7 respectively. The fact that these values are very nearly the same 
confirms the idea that the heavier inert gases may be adequately 
described by a ‘classical’ reduced equation of state (de Boer and 
Blaisse 1948). The present theory, however, gives an expression for 
the magnitude of the constant in terms of other experimentally deter- 
minate quantities. Since no reliable experimental value of y for any of 
these solids is available, we can only check for consistency. ‘Taking 
Robinson’s value for A, c=1-160, we use eqn. (14) to get y=3-56. Using 
simultaneous experimental values for B, 6p, and V,, (Clusius 1936) for 
A at the triple point, the number 3yc(B/62N?/%)°1 can be evaluated ; 
the result is 12-7, as compared with the previously given 12-0 (determined 


or 
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from AH,,/AV,,). The same calculation, using instead the unreliable 
experimental value y=4:5, gives the result 16-2. Such good agreement 
of the two sides of eqn. (20) may be in part confirmation of the theory 
or else fortuitous ; a satisfactory test of the theory in this direction awaits 
accurate y determinations. 
$6 

In the case of solid helium we are in possession of a fairly complete set 
of experimental data with which to check our assumptions, thanks to 
the work of Dugdale and Simon (1953). The specific heat at constant 
volume (,, was measured and the measurements can be represented by a 
universal function f(7/¢), with ¢ a function of volume only. f(x) is not 
exactly a Debye function, however; the equivalent Debye 4’s were 
computed as a function of temperature for various densities. Further, the 
measurements permitted a direct calculation of y by means of the funda- 
mental definition y=—dInd/dIn V; y is constant at approximately 
2-4 over the greater portion of the volume range covered by the experi- 
ments, falling off gradually at lower volumes. The Lindemann relation 
was found to hold within the experimental error over the greater portion 
of the melting curve, not for 4, but for the equivalent Debye 6(V,,, 7’,,). 
If the function ¢(V,,) is substituted in eqn. (2), the constant seems to 
vary steadily along the melting line. Writing (2) in the form 

me2V 2/3 2 


7p == v8 Diz a) eb o0) el towne mire (21) 


where now D is constant, the value n=0-26 is found to represent the 
results of Dugdale and Simon. Thus our implicit assumption as to the 
identity of the two characteristic temperatures occurring in (1) and (2) 
is not fulfilled in the case of helium; the modified Lindemann relation 
(21) must be used for calculations. Further, the following two modifica- 
tions must be made: (a) the zero point energy, which up to now has been 
neglected, must be explicitly taken into account, as it plays a dominant 
part in the behaviour of helium in the condensed phase; and (b) the 
energy expression U—U )=3NkT',, must be altered, since melting takes 
place at relatively low temperatures for helium. 

The zero point energy is included in U,(V), and thus enters the 
Griineisen equation through U,'(V). We can (in theory, at least) separate 
U, in the usual manner, writing U,»=®+U,. Here @ is the static lattice 
energy ; Dugdale and Simon found that good agreement with London’s 
theoretical calculation (London 1936) of ® could be obtained if U, was 
represented by the Debye expression (9/8)Nkdé. The magnitude of d was 
chosen so as to coincide insofar as possible with the equivalent Debye 0 
at 7’=0. Hence eqn. (5) becomes 

ee: ; 9 Nkyd 

Pm=—P Ue 8 am ne vr (U— Uo)r,,, Vin’ stgete (22) 
Now ¢(V,,) is connected with 7, through the modified Lindemann 
eqn. (21) and thus should vary strongly with temperature on the melting 


m 
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line. We tentatively assign the role of the parameter a to the quantity 
@'(V,,), assuming that the melting volumes at low pressure are in the 
vicinity of the inflection point on the D( V) curve. Writing 


GQ = Oe) eames Ee (28) 
eqn. (22) becomes 
9Nk 
Pyta= 5 ee ele iz (Wath, . 2a ee 4) 


The question of the proper theoretical form of U—U, in terms of dhe 
and V,, is a difficult one, for the lattice dynamical theory of solid helium 
is still in an unsettled state. We can expect that U—U, for helium in 
the temperature range covered by these experiments may be adequately 
represented by an expression of the form 


U—U j= 3NkT 2 My. (T/6, . 2... (28) 
j=1 


and as a crude approximation, valid only at the lowest temperatures, we 
neglect all but the first term ; we take M,=74/15 (the Debye expression) 
as giving sufficient account of the correct coefficient. 

We have then, approximately, 


BN ky ho fDiINe, a (Dk\% (m¥2k\3 
Pm o= 8 km =) ae +3 ky 75 Gs hb le 


where 


gi=(1+y)(c— 1), ¢y=(1+y)(1+n)(c—1), 
g=(3y—1(e—1), and ¢=3[1—y(1+n)(e—N]41, . . (27) 


after elimination of V,,. Substituting experimental values for the 
various parameters we find a melting curve equation 


eam DT IL, (9 38) lUe)\ig +e see 2. (28) 


We see that for all but the highest temperatures (when succeeding terms 
give additional contributions) the first term dominates and we get an 
equation of the form of (3), but with c=1-04 as compared with the 
experimental value (Robinson 1952) 1-562. Naturally we should expect 
to be able to fit the melting curve using an equation of the form 


Din pO ET pF KT 6 as) 


for there are more parameters. The fact remains that the melting curve 
of helium is well described (with the possible exception of the extreme 
high pressure portion) by an equation of the form (3). We must conclude 
that the present considerations fail to account for this satisfactorily. 
It may well be that the general procedure set forth here is correct in 
essentials, but that the highly arbitrary step of writing Uo= Or US 
U,=(9/8)Nk¢ is incorrect. A correct microscopic theory of solid helium 
(based perhaps on Born’s recent extension of lattice theory (Born 1951)) 
is at any rate necessary for further understanding of experimental results. 


(26) 
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In conclusion, it is emphasized that no mechanism or description of 
the melting process has been postulated ; indeed, this whole question has 
been by-passed through the assumption of the Lindemann relation from 
the beginning. Rather, these considerations suggest that eqn. (3) might 
be, theoretically, of the same degree of validity as the well known 
eqns. (1) and (2), following as it does quite directly from them under 
certain assumptions. These assumptions seem reasonable in theory ; 
the experimental evidence available for the alkali metals and for helium, 
however, does not wholly confirm them. The evidence is not complete 
for the heavier inert gas solids and predictions have been made from this 
heory which may serve, when further experiments become possible, to 
reject or further verify this explanation of the validity of the Simon 
melting equation. 


ACKNOWLEDGMENTS 


The author is indebted to Prof. F. E. Simon for his interest and 
numerous helpful discussions, to Dr. C. Domb and Mr. T. H. K. Barron 
for advice and assistance in preparing the manuscript, to Dr. J. A. W. 
Huggill for performing the calculations (previously unpublished) of c¢ 
from Bridgman’s data included in table 1, and to the Rhodes Trust for 
financial assistance. 


REFERENCES 


Born, 1951, Festschrift zur Feier des 200 jahrigen Bestehens der Akademie der 
Wissenschaften in Gottingen. 

Ciustus, 1936, Z. Phys. Chem. B, 31, 459. 

CiLusius and WEIGAND, 1940, Z. Phys. Chem. B, 46, 1. 

Davigs, 1952, Phil. Mag., 48, 472. 

De Bokr, 1952, Proc. Roy. Soc. A, 215, 4. 

De Borr and Buarsse#, 1948, Physica, 14, 139, 149, 520. 

Domp, 1951, Phil. Mag., 42, 1316. 

DvuGpaALe and Sion, 1953, Proc. Roy. Soc. A, 218, 291. 

GRUNEISEN, 1926, Handbuch der Physik, 10, 1. 

GUGGENHEIM, 1950, Thermodynamics (Second Ed.), p. 144 (Amsterdam : North 
Holland Publ. Co.). 

LENNARD-JONES and DEVONSHIRE, 1939, Proc. Roy. Soc. A, 170, 464. 

LinpeMAnn, 1910, Phys. Z., 11, 609. 

Lonvon, 1936, Proc. Roy. Soc. A, 153, 576. 

Rosrnson, 1952, Thesis, Oxford. 

Stwon, 1937, Trans. Faraday Soc., 33, 65; 1952 a, Low Temperature Physics 
(London : Pergamon Press), p. 24; 1952 b, L. Farkas Memorial Volume 
(Jerusalem : Research Council of Israel), Special Publication No. 1, p. 37. 

Srvon, RUHEMANN and Epwarps, 1929, Z. Phys. Chem. B, 6, 62; 1930, Z. Phys. 
Chem. B, 6, 331. 

TRUESDELL, 1952, J. Rational Mech. Anal., 1, 125. 


th ye 
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ABSTRACT 


An attempt is made to justify the empirical rules for the spins of odd—-odd 
nuclei recently proposed. It is found possible to do so for certain classes of 
nuclei, which happen to include nearly all those for which the experimental 
evidence is reliable and the rules obeyed. For one of the classes, however, 
it is necessary to make the approximation of assuming nuclear forces of 
infinite range. 


$1. INTRODUCTION 


THE direct experimental measurements of spins of odd—odd nuclei are 
scanty, since they are necessarily limited to stable or long-lived nuclei. 
It is possible, however, to make more or less definite assignménts in many 
cases as a result of studies of decay schemes and ft values. It is naturally 
tempting, in view of the succcesses of the spin-orbit coupling model, to see 
if the values so obtained have any simple interpretation on this basis. 
Empirical rules have been formulated by Nordheim (1951), who summarizes 
his result in these terms : 

“ (i) The individual configuration of neutrons and protons in odd—odd 
nuclei are the same as in odd—A nuclei with the same number of nucleons 
in the odd particle group. 

(ii) If the odd neutron and proton belong to different Schmidt groups 
(i.e. j,=1+4, j,=1+4) their resultant spins will subtract. 

(iii) If the odd neutron and proton belong to the same Schmidt group 
(i.e. 7,=1+4, j,=1+4) their resultant spins will couple to a larger than 
minimum resultant...in many cases near the possible maximum 
Get Jn)e , et. e 

The strict form of rule (iii) ([=j,+),,) will be referred to as (iii’). 
Detailed calculations to be reported later support rule (i) to a limited 
extent. The formulation above is chosen in order to account for the spins 
3 and 4 of the sodium isotopes ?2Na and *4Na, and these as it happens are 
two of the cases in which agreement with rule (i) is least good. However, 
the rule does seem to be generally in approximate agreement with the 
results of detailed calculations. 


* Communicated by the Author. 
+ Now at Department of Atomic Energy, Industrial Group Headquarters, 


Risley, Warrington, Lancs. 
2C2 
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Most of the possible states of an odd—odd nucleus have spins of the order 
(jp t+jn), and in the absence of any special symmetry requirements such 
as those applying in even-even nuclei, we would expect spins of this 
order. Rule (ii) and rule (iii’), therefore, are more significant than rule 
(iii). Other things being equal, and accepting rule (i), it is clear that 
spin-dependent forces will favour parallel spins of proton and neutron, as 
in the deuteron, and so favour both (ii) and (iii’). | But nuclear forces are 
predominantly non-spin-dependent, and it seems worthwhile to seek 
theoretical justification of the rules more searching than the remark just 
made. This is attempted in the following ; but we must first point out 
that the experimental evidence for the rules is in fact more limited than is 
claimed by the proposer. 

Nordheim (1951), indeed, quotes a long list of spins, with shell assign- 
ments for the odd particles. Many of the latter seem a little arbitrary, 
and indeed, for A >60, there is, in general, sufficient competition between 
shells for almost any datum of the type available to be fitted into the 
scheme of the rules. We thus need only consider nuclei with A<60. 

Now nearly all nuclei with A<60 fall into one of the classes detailed 
below, for which we shall justify the rule theoretically, and table 1 shows 
that those which do not, frequently do not satisfy the rule also; especially 
when it is remembered that rule (i) has been tailored to fit the Na isotopes. 


Table 1. Odd-—Odd Nuclei, A<60, Showing Agreement with Rules (ii) 
and (iii’) 


Class A Class B Others 


Rule obeyed or 2eR tN a 16N, 22Na, 24Na 
he 50 


| Rule possibly eC Leo aro 0p. S8P py tos hee ee Ned 
| obeyed 


Rule not obeyed | ®Li,* 3®C] F: = 26A], 40K 448¢, 
46S, 2Mn, 48V§ 


Insufficient 
»vidence = ts 20 2 
evidence F, 28A] 


* SLi is not a true exception. In this light nucleus the tensor forces must be 
considered. See Hitchcock (1952). 


§ See Nussbaum, van Lieshout and Wapstra (1953). 


§2. DETAILED CALCULATIONS 
The two classes of nuclei mentioned above are : 
A. 1 nucleon or 1 hole in both proton shell and neutron shell, protons 
and neutrons in the same shell. 
B. 1 of the odd nucleons is in an s,/) shell, and there is 1 nucleon or 1 hole 
in the shell of nucleons of the other kind. 
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In Class A we are only able to justify the rules on the assumption that. 
the range of the nuclear forces is infinite, and consists of Majorana, Wigner, 
Heisenberg, and Bartlett forces only, the Majorana and Wigner forces. 
dominating. 

The approximation of forces of infinite range is, of course, poor, but in 
the absence of anything else, these calculations may not be worthless. 
especially as the calculations reported in § 3 show that the result continues. 
to be true for finite range forces in some special cases. 

We shall denote by (a) the wave-function of a single particle, magnetic 
quantum number a, in the shell nl, ; the wave-function of a hole, which is a 
Slater determinant of order 2j, will be denoted by (a’). The wave- 
function of the state of spin J (with M=T) of two unlike particles (or 
particle and hole, etc.) is a sum of products of the functions (a) (or (a’)) and 
will be denoted by 


1 
PL) ay (Jk j+)} 
2;—1 
(or ae a, {((j—k)(I—j+k’)}, etc.) 
(The «, are Wigner coefficients; we do not require their exact. 
form.) 
AG NE 
k 


(The closed shells inside the nucleus can, of course, be neglected.) 
We shall denote the energy integral 


$ (a),* (b)o* Plc)y (d)_ 27, Bra, 


(where P is an exchange operator ; W, M, H, or B, in the usual notation) 
by (ab) (cd)p. We shall consider one case in detail, and quote the results. 
for the others. 
Class A(1). j=Il+4, one particle of each sort in the shell. By direct. 
calculation : 
(ab) (cd) 7=8 ae Soa , | 
442(ab) (cd) p= Sac 8a (2)* +206) 
Si Omae 1 On) dat {G+a)(j—-4+ 1)G—8)G-+6+ AW Ets (1) 
#8) oop, air ee HT 8) D= A), 
(ab) (cd) y=(4b) (de)z, 
(ab) (cd)y=(ab) (de). | ; 
I, the total spin of the nucleus, can take every value from 0 to 2), and since 
(Qj —I—b)(E-+k-+ D)}? at (RAV QG—B)}Y ay=0 © (2) 
so that “,=(—1)*1 Hoi T—k Te ws oe) 
Thus the binding energy Hp(/), of the state of spin J due to the forces P is 
in arbitrary units 
2,—12,—I 


Ep(I)= 2 2 am, (J—k I—j+k) (G—m 1—jrm)p. 


k=0 m=0 
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Hence F,, (/)=1, from (1). 


2,-1 


4jBy (I= 2 oy [tas-r-v {(2—HL—k) +k (23 -1—8)} 
h=0 
foes -rner {e(+)(2J7—k+1)(23—-I+k+ 1)? 
+951 h-1 EN ae aa 

—(—1)'*! (7+ 1)—2)), from (2) and (3). 
47°E , (I)=(1(f+1)—29), 
Ey (1)=(—1)™*. 
Since the Majorana forces dominate, the ground state is that for which 
E,, (J) is the largest, i.e. /=27, in accordance with (iii’). 


Class A(2). If, instead of two particles in the previous case, we have 
two holes, the results are exactly the same. 


Class A(3). If, however, we have one hole and one nucleon, we obtain : 
Ey (I)=—1, Fe 
Ey L)=CU+1)—2))/49?, 
Ey (1)=0, 140, 1: =a complicated expression <0, J=0, 1. 


Ey (1)=0, 140: =—2, I=0. 
In this case only the Bartlett forces separate the states of high spin, and 
these favour rule (iii’), provided the coefficient is positive, which is so for 
all proposed forms of the total P. 

Class A(4), (5), (6). Ifj=l—3, similar results hold. 

In Class A, therefore, the rule is justified in the strict form for the 
infinite range case. Since we are concerned only with light nuclei, this is a 
not impossible approximation. 

For Class B nuclei, the demonstration can be made for a general radial 


dependence of the nuclear forces. We shall suppose that it is the second 
particle which is in a state of spin}. Then 


#(5+3)=()i Be 

+ (3-4) =[G—Vx (B)2— V29(9) a (—H) ol /(29++1)?, 

Hp (J+2)=(2) (9 

Bp j—qy— PLE MG=1 Y= 2V BGR NG 4G — DUG =D 
(2)+-1) 

Now, imagining the indicated operations to be performed under the sign of 


integration, and remembering that the conjugate complex of the first two 
wave-functions is taken, we have :— 


GO“ D=G DEAN De U Rai 2} 
aie 2)(9 2)p—( “i Oe: Des V2) 
and similarly, 
(J—1 3)G—1 2)p=[(2—1(5 39 Het (5 —4)(5 —4)p)/2i 
so that AB» =Ep(j-+4)—Ep(j—H= (GG p—G GF 
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Then, writing down the integrals, 
4H y,—=AE,=0 in all cases; and we have also :-— 
Class B(1). j=1-+4. 


2j+1 
AE,= j : | ¥4(0,) ie |B (73) ? [R.’ (72) I? V(142) dr, arg, 


» 
~ 


+1 
Nie } Yi8() ¥'(65) Ryg(?y) Rygls) Ryy ol) Ry alt) 


2) 
x V (14) Br, Br, 


In the integral for 4H,, the integrand is positive. Hence 4H,,>0. 
Also 


27-1 F 
4E,+ 4E y= a f | 21s) Yi(0;) Ry (TatRyrjl(ro) Y'(0)* Ry (11) ; 


4 
KVetas) Creer, 
which is positive for the same reason. Hence H,>|E,|. The resulting 
Hp is thus certainly positive of P has either the form proposed by 
Rosenfeld (1948) (=0-93M —0-13W -+-0-46B—0-26H) or that proposed by 
Serber (P=0-4M-+0-4W-+-0-1B+0-1H). Thus the ground state certainly 
has J[=j-+4 in accordance with the rule. 

This is the simplest case: we shall not write out the integrals for the 
other cases, but express the results in words. 

Class B(2). j=l—4. 

By an argument more complex than that above, but still elementary, 
we can now show that 4H#,<0. Also —4H,>|4EH,|*, so that, as before, 
we have [=)—4, in accordance with the rules, for all forms of V and P. 

It is to be noted that the proof is possible in case B, only because of the 
failure of the dominant Majorana and Wigner forces to separate the two 
states, so that spin-dependent forces, as in the deuteron, favour parallel 
spins. 

$3. Errects oF FintrE RANGE oF NUCLEAR FORCES 

For forces of finite range, a general calculation in case A is no longer 
possible. Calculations have been carried out for the nuclei *F and ”F on 
the (false) assumption that these nuclei have orbitals in the d;,, shell only. 
The technique is described in a following paper (Hitchcock 1954). It is 
found that, except for very short range forces (iii’) continues to hold, but 
that for these the calculated spin of F is 0, and of #*F, 4. Similar 
calculations on Sc*? and Sc*® together with those of Kurath (1951), and 
those of Flowers and Edmonds (1952) suggest the following rule :— 

For cases A(1), (2), (4) and (5), (iii’) holds, except for very short range 
forces, or large nuclear radii; for cases A(3) and (6) rule (iti’) holds for 
long range forces (or small nuclear radii), but for forces of shorter range (or 
larger nuclear radii) the calculated spin is (27—1), so that (iii) holds. No 
general proof of this proposition, however, is available. 


ie eee ao SY et ee ee ee 
* This inequality cannot be proved with the same generality as the others, 
but is certainly true for all the usually assumed forms of V. 
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For forces of sufficiently short range, Pryce (1952) has shown that, in 
shells of large j, the least possible spin is probably always preferred. All 
the calculations mentioned above, however, indicate that this ceases to be 
true at the range of forces which is of interest, at any rate in the light 
nuclei considered here. 

An attempt to extend Class B to include nuclei containing a py/> 
nucleon fails, since we no longer have 4H#,,=0. A similar proof to those 
of class B does, however, readily apply to the two extreme cases of infinite 
range and zero range forces, since for these we do have 4H ,,=0. 

We conclude that the agreement between the proposed empirical rules. 
for the spins of odd—odd nuclei and experiment is partially accounted for. 
We have already pointed out (§ 1) that the spin-dependent forces must 
always favour the truth of rules (ii) and (iii’). If there is indeed a tendency 
for the rules to be obeyed in other classes of nuclei, this may suffice to 
account for it. Certain attempts to extend these general calculations 
have proved unsuccessful, and more detailed calculations (Hitchcock 
1954) indicate that the rules no longer hold theoretically. 


The author’s thanks are due to Professor M. G. Mayer, and to 
Mr. J. M. C. Scott, for valuable criticism and discussion. 
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XLVI. Spins of Odd-Odd Nuclei in the j-) Coupling Model: IT 
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[Revised MS. received December 14, 1953} 


ABSTRACT 


The spins of a number of odd odd nuclei are calculated on the basis of 
the spin-orbit coupling model. The study is limited to light and medium 
nuclei for which experimental data are available. The agreement between 
theory and experiment is generally good. No systematic rules are pro- 
posed, and it appears that no simple ones exist, except in the cases 
previously treated. 


§1. IntrRopuUCcTION 


It has already been shown (Hitchcock 1954) hereafter referred to as I, 
that the empirical rules proposed by Nordheim (1951) for the spins of 
odd—odd nuclei can be justified in certain restricted cases. It is the object 
of this paper to extend the investigation of such spins to other nuclei, 
in order to see if the rules may be more generally applicable. (The 
rules are quoted in I, and will hereafter be referred to as (i), (ii), (iii) 
and (ii1’).) 

The j—) coupling model is taken as a basis for all the work, and, except in 
one case, it is assumed that nuclear forces are entirely two-body in char- 
acter, though different exchange operators are considered. We also assume 
charge independence of nuclear forces, that is, that proton—proton, 
proton-neutron, and neutron—neutron forces all have the same form; and 
neglect the Coulomb force entirely. (It has been found by trial in a few 
cases including °°V, one of the heavier nuclei studied here, that its inclusion 
does not affect the results provided that no variation in the core configur- 
ation is allowed.) 

We shall refer to the closed subshells of a nucleus as the core, and the 
remaining nucleons will be called orbital nucleons. In all calculations of 
this type we are concerned with the lowest excited states only. Con- 
figurational mixing, both in the core and in the orbital nucleons is neglected 
and the core itself can therefore be neglected, except in so far as it 
provides the source of the fictitious potential well in which the orbital 
nucleons move: for on account of the complete spherical symmetry of 
closed shells, the energies of orbital nucleons in the field of the core are 
equal, provided no orbital nucleon is excited into another shell. Thus we 
need only consider the interactions between the orbital nucleons. 
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From this, it follows without difficulty, that the low-energy level schemes 
of what we here call conjugate nuclei are identical, apart from small 
effects due to differences in nuclear radius. The nuclei conjugate to a 
nucleus with NV, protons in the n4/,j, shell and N, neutrons in the 741), 
shell are defined to be the nuclei with (2j,+-/— N,) protons (or neutrons) 
in the n,/,j, shell and (2j,+-/— NV.) neutrons (or protons) in the 7/57. shell. 

In many cases, the level schemes of conjugate nuclei are dissimilar. 
This must be ascribed in part to the effects of different nuclear radius, or 
to the Coulomb force, but these are too small to account fully for the effects, 
which are therefore due to partial failure of the j-j coupling model itself. 
(We shall say that the model ‘ fails’ if the effects of configurational 
mixing become large enough to affect the results.) 

In general, indeed, these calculations, and the comparison of their 
predictions with experiment, provide us with a sensitive test of the extent 
of the validity of the assumption of spin-orbit coupling. In an odd-odd 
nucleus, which has none of the simple symmetry properties of other kinds, 
there are many possible ground states, between which we cannot choose 
without explicit calculation. It is thus satisfactory that only one 
definite disagreement with experiment is found. 


§ 2. ForcES AND RapiAL WAVE-FUNCTIONS 


The greatest uncertainty in all calculations of this kind arises from 
our lack of knowledge of the form of the two-body forces between 
nucleons. For want of better information we have assumed that the 
inter-nucleon potential is of one of the forms : 

(i) V(r)=V, exp(—ar?)P, referred to as G 
succeeding, where 


, for Gaussian, in the 


V,=42 Mev, 
a%=0-1089 x 1026 cem-?. 
(ii) V(r)= Vo (exp(—zr)/r)P, referred to as Y, for Yukawa, in all the 
succeeding, where 
Vo>=89 x 10-18 Mev cm, 
“a=0-858 x 1038 em7!. 
(iii) V(r)=V,3(r)P, referred to as 8 in all the succeeding. 
The numerical values are taken from those derived by Blatt and 
Jackson (1950). P is supposed to be an exchange operator: we take 
(a) P=0-:8M-+0-2H referred to as MH, or 
(6) P=0-4M-+0-4W+-0-1B+-0-1H, referred to as S, for Serber, in the : 
succeeding. 
The two forms are indistinguishable in case 8. 
The tensor force has been neglected. In a large nucleus, where its 
range is comparatively small, this approximation may be justifiable, 


but as the calculation on Li (¢.v.) shows, it would be desirable to consider 
it. It is hoped to do this in future work. 
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Further uncertainty arises in the form to be adopted for the radial 
wave-functions. In theory, a sufficiently accurate approximation to 
these may be obtained by a variational calculation, but in practice it is 
also necessary that they be in a form suitable for the performance of the 
integrations arising. We have therefore chosen in case G 

R,.i(r)=exp(—Br?)r"H,,,,(Br?)—* harmonic oscillator’ functions and in 
case Y 

f,,,(7)=exp(—fr)r'L,, (Br) — hydrogen ’ functions, where in each case 
the H and L are the suitable polynomials. 

f is determined, following Kurath (1951), by fitting the wave-function 
as well as possible to the corresponding wave-function of a square 
well of radius 77>=1-48x10-1341/3 em. In case § the form chosen is 
irrelevant, provided configurations with nucleons in varying shells are 
not considered, as we shall not. Since this method is uncertain, a range 
of values of 6 has been used in all the calculations, and it is found that 
in most cases the order of levels is not sensitive to the value chosen. 
However, the actual level separations depend very sensitively on B, 
so that it is not possible to predict the values of the separations to an 
accuracy sufficient to make comparison with experiment worth while. 
Where the calculated spins differ for larger and smaller admissible 
values of the nuclear radius, we shall distinguish the case of the larger 
radius by a suffix 1, and the smaller by a suffix 2, thus: G,,YWH,. 
If one of these is less likely than the other it will be bracketed. 

We may note that the wave-functions of case Y are probably less like 
the actual ones than those of case G, and the use of the former, when two 
or more shells are involved, leads to violent anomalies (Hitchcock 1953). 
Case G is therefore to be preferred if the results differ. In fact, as we shall 
see, case GS nearly always agrees with experiment. 


§ 3. CLASSIFICATION OF Many-PartTIcLE WAvE-FUNCTIONS 


It may happen, and in general does, that a number of orbital nucleons 
can have more than one configuration of the same spin and isotopic spin. 
A further classifying observable is, therefore, required. Ideally it is 
desirable to find a third observable which is actually conserved under the 
inter-nucleon forces, and which, like spin and isotopic spin, lends itself 
readily to group-theoretical treatment ; but no such observable seems to 
exist. It is necessary therefore to consider a linear combination of the 
wave-functions arising in the classification chosen, calculate the energy, 
and minimize. A classification which approximately satisfies the ideal 
condition reveals itself by small cross-terms in the parametric expression 
for the energy. 

Flowers (1952) has introduced the concept of seniority, which satisfies 
the required conditions almost exactly. This is certainly the most 
elegant, and probably the best classification proposed. Kurath (1951) 
obtained wave-functions which contained a minimal number of terms, 
subject to being eigenfunctions of spin and isotopic spin. For the simple 
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cases he treated this is adequate, but the method is of little use in more 


complicated cases. ; 
As the work of Flowers was unpublished at the time these calculations. 


were begun it was decided to adopt a classification in which proton spin 
and neutron spin were separately conserved. This has the following 
advantages : 

(a) It also approximately satisfies the ideal condition. 


(b) It lends itself to physical interpretation, and enables us to see at 
once to what extent rule (i) is valid. 

(c) It breaks a complicated calculation into a number of small sections, 
making checking easier, and reducing the likelihood of computational 
error. 

(d) It is found by trial that states of very large proton and neutron 
spins occur to a negligible extent in the low excited states. This enables 
the neglect of these states, and consequent great simplification of the 
calculations. 

(ec) It should be noted however that the states obtained are not eigen- 
states of the isotopic spin. 


§ 4. SprciFic NUCLEI 


We shall consider conjugate nuclei together, where both exist, and 
summarize the results only, without entering into the details of the 
calculations, or quoting numerical results. We shall not treat nuclei in 
which one nucleon is in an s1/, state, for which a general proof has been 
given in I, unless some special consideration makes it desirable. The 
quoted experimental results are derived from Way et al. (1950). 


SLi and 1°B 
Measured spin of ®Li 1 
Measured spin of 1°B 3 
Calculated spin of ®Li l 
Calculated spin of 1°B 3 


The work on *Li has already been reported (Hitchcock 1952 a). In 
this case alone the tensor force was considered. It is shown that it is 
possible to account for the spin 1 by assuming strict j-j coupling, although 
the contrary has been asserted. The effect of the tensor force, however, 
decreases rapidly with increase in nuclear radius, and a rough calculation 
for 1°B indicates that in this nucleus its effect is not sufficient to invert 
the levels of spins 3 and 1. The position with regard to excited states is 
equally satisfactory. ®Li has two excited states of spins 3 and 0, which 
occur in the order predicted. The two lowest excited states of 1°B have 
spins 1 and 0 according to both theory and experiment. 

It should be pointed out that in spite of these satisfactory results, 
it is doubtful if there is, in fact, true j-j coupling in ®Li. 
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Much work has been done on these nuclei assuming intermediate 
coupling but neglecting the tensor force (Inglis 1953, Racah and Zeldes 
1951, Zeldes 1953). It seems probable that the effects of both inter- 
mediate coupling and of the tensor forces should properly be taken into 
account. The work of Zeldes seems to suggest that sufficiently large level 
Separations cannot be obtained on the basis of intermediate coupling 
alone. (Inglis’ calculation assigns the nuclear parameters to fit the level 
schemes in a manner which seems to allow different values of these 
parameters in the two nuclei.) 


14 
Measured spin i) 
Calculated spin 1 


This nucleus was considered, allowing for the possibility of p,/>—Pps/o 
mixing, in order to see if the remarkably long f-half-life of 4C, and the 
associated lack of the ground-state transition in 0, could be accounted 
for by ‘ accidental’ vanishing of the nuclear matrix element. A tensor 
force, of the same range as the ordinary forces, was assumed. In fact 
too much mixing is found, and the ft value is increased in consequence, 
assuming Gamow-Teller selection rules, by a factor 3 only. 

The predicted magnetic moment is 0-63. If ps. mixing is neglected 
we get 0:37. The measured value is 0-40. This suggests that if the 
mixing were less for some reason (e.g. if the ps/.—Ppy/, Separation is 
greater in 4N than in “Li, the value chosen) there would be better 
agreement. But in fact no improvement is found if the amount of mixing 
is adjusted to fit the magnetic moment, though consideration of the 
mixing of higher states could improve the fit. 


18s 
Experimental spin tror.0% 
Calculated spin, Y and G, i 
Calculated spin, G, and 6 0 


Here the possibility of 1d;/.-2s,/, mixing was considered. As already 
reported (Hitchcock 1953) the pairing energy of the 2s,/. nucleons is 
greater than that of the 1d;,, nucleons, for cases of long range forces. 
Consequently, the ground and first excited (calculated spin 0) states are 
predominantly (s1/2)? configurations, the lowest (d;,/,)? configurations 
coming much higher. This may account for the lack of low excited 
states in this nucleus (Hitchcock 1953). 

2Na 


Measured spin 

Calculated spin, 5, G and (YS;,) 
Calculated spin, Y MH 
Calculated spin, YS, 


w ord w 


* In view of the 0->0 transition in 140, a not dissimilar nucleus, the value 0 
cannot be excluded. 
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Agreement with experiment is unsatisfactory. 281/, mixing has been 
neglected, but a rough calculation indicates that its inclusion is unlikely 
to improve the situation. The f-transition to an excited state of *Ne 
has a log ft value of 7. Since the shell model does not allow for parity 
change here, this is exceptionally large. The calculations do not uphold 
the suggestion that this is due to the re-arrangement of the configuration 
(ds/)3-3/a, a8 rule (i) suggests, since 30% of the wave-function corres- 
ponds to the configuration (d5,2)—5/2(d5/2)}—5/2: 

It is possible, however, though unlikely, that s,,. mixing occurs to such 
an extent that not only is the spin 3 state depressed to the ground state 
but that the B-transition is ‘/-forbidden’ (vide *P). This solution is 
advanced with little conviction. 


244Na and 29 


Measured spin of 74Na 4 
Experimental spin of 27°F 2 or 3 (22) 
Calculated spin, GMH, YS, YMH, 5 
Calculated spin, 6, GS, (YM#H,) 4 


2814/2 mixing has again been neglected. The experimental indications 
for 2°F are conflicting, so that disagreement with experiment is not 
established here. For *4Na agreement is excellent. The level schemes of 
these nuclei have been extensively investigated. There is no evidence of 
the correspondence expected in conjugate nuclei. 

281/, mixture in 2°F, 2Na, and *4Na 

The results for SF suggest that the approximation of neglect of 2s,/, 
mixture in the preceeding nuclei may be very poor. It is of interest to 
proceed to the opposite extreme, and suppose—as indeed the calculation 
for 18F suggests—that there is complete replacement of d;,, pairs by s4/> 
pairs in these cases. Calculated results now become : 


°F, Configuration ((81/.)?d5/2), dsj. 


Calculated spin, 5 0 

Calculated spin, G and Y 5 
Configuration ((81/9)ds5/s), 84/9. 

Calculated spin 3 (see I). 


#2Na, Configuration d;/5, d5/>. 
As 2°F, first above. 


**Na: The new configuration (d5/9)3(d5/2) is the conjugate of the former 
(d5/2)°(d5/2)*, so that the results are unaltered. In this case, therefore, 
intermediate mixture will also probably lead to the same result. 

There is little improvement ; and in any case, the spin 3/2 of 24Ne shows 
that its configuration is not ((8y/9)? 5/2), (84/9)? as the assumption above: 
would require. 

32p 

Id5/9~281/2 mixing was considered. The log ft value of the B-transition 

to the closed-shell nucleus *8 is 7-9, and the Kurie plot is straight. The 
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shell assignments do not allow of parity change, so this is presumably an 
allowed transition. Nordheim (1951) classifies it as ‘J-forbidden’; 
suggesting, in effect, that transitions with 4J=1, Al—2 (S1/2—As9, ete.) 
cannot be placed in the same class as other allowed transitions ; and indeed 
the Gamoy—Teller allowed matrix element vanishes for such a transition. 
The observed decay must, therefore, be due to a ds. contamination in 2P, 
which must be of the order 0-1% to account for the observed half-life. 
For a suitable and admissible choice of the nuclear radius this can be done 
in both cases G and Y, and at the same time a separation between ground 
and first excited states of the right (70 kev) order of magnitude obtained. 

It is pertinent to ask why a precisely similar state of affairs does not 
apply in **P, which has a log ft value of 5. The answer lies in the fact that 
the interaction between the (d3;.)~! neutrons and the (d;,;,)~! protons in 
the ground (spin 1) state of **P is much greater than that between the 
single d3,, neutron and the (d;,/,)~1 protons in 32P, whereas the interactions 
of the two neutron configurations with an s,/, proton are the same. Thus 
the d;,, contamination in *4P is appreciably greater than in ®2P, and the ft 
value of the f-transition correspondingly less. Thus the ad hoc assign- 
ment by Nordheim of the fifteenth proton to an s,/, state in 9?P and toa 
d;,. state in 34P has approximate justification. 

Another explanation is possible. If the spin of **P is in fact 2, as is 
indicated theoretically as a possibility, then the proton configuration is 
mainly (d;/,)-1 and we have a second-forbidden transition with A/J—2, 
Al=0. Following Nordheim’s nomenclature we may call such transitions 
‘/-preferred ’, since calculation of the matrix element indicates that 
the ft values of such transitions should be reduced by a factor of the 
order 3 < 10-3, which leads us to expect log ft values between 7-5 and 8-5 
for such transitions. This fits the value for **P, but the allowed shape 
of the B-spectrum makes this explanation less plausible than the previous 
one. 


36C] 
Measured spin 2 (3 2) 
Calculated spin, GS, 6 2 
Calculated spin, Y, GMH 3 

Agreement is satisfactory. 

38] 
Experimental spin Os ror? 
Calculated spin, 6 2 


The shell model allows no spin less than 2, so it is satisfactory that 2 is 
found. It is the most tightly bound by so much that it is likely to be so in 
cases @ and Y also for all reasonable values of the nuclear radius. The 
first excited state is predicted to have spin 5, which suggests isomerism, 
but a finite range calculation, will be required to establish if the model 


really predicts this. 
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40K 
Measured spin 4 
Calculated spin, 6 4 
According to Kurath (1952) the calculated spin is 4 in case G also. 
44Sc and **Mn 


The calculations (case 5 only) indicate that the spins of the low excited 
states, in ascending order, are | (ground), 6, 4, 5.* 

These nuclei both show isomerism and, according to Goldhaber and 
Sunyar (1952), the isomeric transitions are M4 in **Mn, and M4 or E4 
in “Sc ; an E4 transition implies parity change, which is not acceptable 
to the shell model. M4 is therefore assigned to both transitions. This 
is slight confirmation of the theory of conjugate nuclei. 


Decay scheme of “Sc . Assigned Decay scheme of *?Mn. 
values of spin and parity given 
in brackets; (all.) indicates 
allowed transition. 


That the ground state is predicted to be of low spin is satisfactory, at 
first sight, and we must assign spin 5 to the second isomer. The 
depression of the state of spin 5 below the (very close) states of spin 4 
and 6 can very probably be due to the effect of finite range nuclear 
forces—a similar depression is predicted in *8Se (q.v.). We therefore 
assign spin 1 to 4-hour “Se and 21-minute °*Mn*; and assign spin 5 
to 2}-day “Sc* and 6-day °*Mn. It is difficult to see, if these assignments 
are correct, why the f-transition to the ground state of the daughter 
nucleus is not observed, but the lowest state of spin 2 is predicted to be the 
seventh excited state, and cannot be depressed to the ground state by 
considerations of finite range. Agreement with experiment is thus not 
satisfactory. 


* Dr. Flowers informs me that there is an error in the calculations reported 
in his paper (1951) which disagrees with those reported here. 
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46Sc and °°V 


Measured spin of 5°V 6 
Experimental spin of 46Se 4 or 5 
Calculated spin, 8 6 
Calculated spin, G and Y 5 


The work of ®°V has been reported previously (Hitchcock 1952 b), and 
the same calculation applies for 4*Sc. It is certain, that for all nuclear 
radii in the relevant interval the calculated spin will be either 6 or 5, 
and it may be that the slightly larger nuclear radius of 5°V is adequate 
to ensure its conformity with case 5, whereas it is necessary to invoke 
finite range forces to explain the spin of 4*Sc. The difference in radius 
must be small, however, and the calculations do not indicate that the 
cross-over point is in the region of interest. 

*6Sc shows E3 isomerism (Goldhaber and Sunyar 1951). The calcula- 
tions suggest spin 2 for the isomer, which is in accordance with the decay 
scheme. No such isomerism, however, is apparent in ®°V. 


§ 5. CONCLUSIONS 


We conclude that the assumption of strict jj coupling is not inconsistent 
with the observed spins of light and medium nuclei, especially, it 
will be noted, in case GS. There is no reason to suppose that it will not 
also account for those of heavy nuclei, but there are few experimental 
data for comparison, and shell assignments are uncertain, so the laborious 
calculations are not profitable. 

Agreement with experiment is generally satisfactory (except for ?2Na, 
and to a lesser extent for “Sc and °*Mn) so far as spins are concerned. 
The agreement of the level schemes of conjugate nuclei is not good, 
except for the 4/=4 isomeric transitions in “Sc and °**Mn. A tendency 
for the spins to increase through a shell, pointed out by King and Peaslee 
(1953) is, of course, contrary to the j—j coupling model, and could not be 
justified by these calculations. 

There is no theoretical justification of rules (ii), (iii) or (iii’). except in 
the cases considered in I. 
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ABSTRACT 


If an insulating liquid is placed between the plates of a condenser the 
electric field gives rise to induced intermolecular forces which are 
proportional to the square of the field strength. On account of these 
forces there is a change in the radial distribution function. The field 
effect on the distribution function is calculated for the case of thermal 
equilibrium and for steady flow. The corresponding increment in 
viscosity compares not unfavourably with experimental facts. 


$1. IyTRODUCTION 


THERE is, so far, no molecular theory of the effect of an electric field 
on the viscosity of liquids as observed by Andrade and Dodd (1950). | 
In this paper an interpretation of this effect is attempted by considering 
the intermolecular forces which are induced by an electric field and their 
contribution to the flow of momentum. Ignoring some finer aspects 
a simplified molecular model is used by assuming that the molecules 
interact with central forces on which the forces of a molecular dipole 
may be superposed. The interaction of the molecules in an electric field 
is then obtained without any further assumption. The contribution of 
the induced forces to the coefficient of viscosity is arrived at by applying 
the theory of transport processes as developed by Kirkwood (1946) and 
Kisenschitz (1949). 


$2. INTERACTION OF MOLECULES IN AN ELEectrrRic FreLp 
The average electric moment of a molecule is equal to 
M=[a+(w/3kT)|F=oF +My... .. . (I) 
(Debye 1912) where « and y are the electronic polarizability and the 
permanent moment; F is the field vector. 
Molecules in an electric field should accordingly interact like parallel 
dipoles with a potential energy of interaction 
PCE) ae cee en OL 
where K(r)=K(r, 8)=(1—3 cos? $)r-2 |g wt (3) 


r being the distance vector. 


* Now at the Admiralty Research Laboratory, Teddington. 
+ Communicated by the Authors. 
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Ignoring at first the first term in (1) the validity of eqn. (2 a) is readily 
established. The potential energy of a pair of polar molecules in an 
electric field is equal to 


Xr=Xoo(”) +[7*(uy » Ye) —3(uy . )(We . Fr) >— (uy +p) . F 7 (4) 
where 4, U2 (uy=e—=) are the molecular dipoles ; yoo is the familiar 
potential energy of molecular interaction. An exact knowledge of the 
potential energy at small distances—where eqn. (4) does not apply—will 
not be required in this paper. 

The force of interaction is —grad y,, where 


xXu—=—kT log [(1/167%) | exp (—xp/kT) d(pylu) U(we/u)] . . (5) 
By expanding to the third power of the exponent it is found that 
Xu=Xoo— (wa) 8kTr®) — (u2F?/3kT!) + M PK (r, 9) 
which is in accordance with (2a). 

The interaction of non-polar molecules is calculated by applying 
quantum mechanics to the electrons. The second order perturbation 
of energy is familiar and equal to —a«F?+yo9. In the third order an 
expression is found which is very nearly equal to («F)?K(r.%) in 
accordance with eqn. (24). 

The complete interaction of polar molecules is—for the simplified 
molecular model—obtained by successive application of perturbation 
theory to the movement of electrons and evaluation of the angular 
partition function of the molecular dipoles. The result is in agreement 
with (2a) and conveniently written in the form 

VareNor® far Kh (re) as. | (28) 
where yy is the potential energy in absence of a field and «’ is the mean 
polarizability. 

In this paper eqn. (2 5) is applied to a pair of molecules in the interior 
of a liquid. J is then the internal field strength and related to the 
external field strength in terms of the dielectric constant « according to 


FEW) (cee) eee en! (6) 


$3. Liqgurp In THERMAL EQUILIBRIUM 
If fo(r), f(r) are the distribution functions in relative coordinates in 
absence or presence of an electric field, ‘ effective * potential energies of 
interaction ®), ®y are defined by writing 
f=exp (— Ok) 
and the negative gradient of the effective potential energy may be 
regarded as the force determining the relative movement of the molecules. 
The effect of an electric field on the properties of a liquid is determined 
by its effect on ®p. | 
A rough estimate of this effect is made by assuming that the difference 
of ®,, and xy+®, arises from the interaction of a shell of parallel dipoles 
with a molecule when another molecule is located at the centre of the 


2D2 
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shell. Let the radial distribution of the dipoles in the shell be specified 
in terms of a distribution function f,(7) then the effective potential is 
found to be 


Dy = Dy xa} Se MRO) {3 [Psto'ye'® dr’ | fale ar’. Sara) 
0 r 


Using a distribution function f, of simple analytical form and 
corresponding to 10 nearest neighbours, expression (7) was evaluated. 
The result was applied in numerical tests, as specified in §6, but not for 
calculating quantities to be compared with experiment. 

Improving on this calculation, consider the energy of interaction of 
the molecules labelled as 7 and j 


Xults J)=Xal 3) =Xol" es) +ICK j; 
where K,, is written for K(r,;). The probability distribution for the 
simultaneous positions of the molecules 1 and 2 is given by the expression 


full, 2)=(0/2) | exp (QUT) Exalij)) drys. ay 


49 
where V is the volume, Z(7', V, ) is a normalization factor satisfying 
the relation [M~1(0Z/0M)|y,_»>=0; in the integrand molecule 2 is 
placed at the origin of the coordinates. 
Expanding f,, in powers of M?, the coefficient of M? is equal to the 
required increment of the distribution function. It is equal to 


KI'(bu—$o) fol” 2=K oxfoll'a)+2N | fo (1, 2, 3)K 93 digs 
+N? | fy (1, 2, 3, 4)Kg4 digs dry. = a om 

Here f, and f)? are the distribution functions for triplets and 
quadruplets of molecules and N is the number of molecules present. 

Equation (8) does not depend on any approximation other than implied 
in the use of a constant internal field. The right-hand side of (8) can, 
however, be evaluated only if the unknown distributions are eliminated 
conveniently by means of the ‘ superposition approximation ’ 


Jon (1, 2, 3)=fo(1, 2) fo(1, 3) fo(2; 3), 

fore (1, 2, 3, 4)=f, (1, 2, 3) fo(1, 4) fo(2, 4) fo(3, 4). 
By entering with (9) into (8), ®), is expressed in terms of known 
distribution functions. 

The last term in (8) is, in spite of the factor N2, much smaller than 
the preceding term on account of the angular function in the integrand 
which reduces the integral to zero unless the four molecules are situated 
close to each other. It was accordingly possible to neglect this term. 
The remaining terms are reduced to the following expression 


Dy =P,—2M*P,(cos 9)[r2+2n7NJ(r)]  . 2 ©. (10) 
where J(r)= | | Solr’) fo(r’)r’—1 sin 6 dé dr’ 
0 Jo 


and 7’ =(r?-+-r'2— 2rr’ cos 6) 12 


(9) 
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and P, is the Legendre polynomial. Expression (10) may be regarded 
as a better approximation than (7) but its evaluation requires extensive 
numerical work. Expressions (7) and (10) can be written in a common 
form in terms of two functions of 7, denoted by @, and @, 


Py=O,+ M7 O,+26,P(cos9)] . . . . . (11) 


which will be used in the following calculations. 


§4. Tue Sreapy Non-Untrorm State 


The distribution f4(r) which is appropriate to a liquid in steady flow, 
is determined by the ee equation 


Of 4[Ot—=div {(kT ym) grad f4+[(grad Jyy/ym)—e]f4}—=0. . (12) 
Here ¢(r) is the macroscopic velocity, m is the molecular mass and y a 
friction constant. Let 


Cre=0, ==) Cy=Az 


where A/2 is the constant yz-component of the rate of shear. Neglecting 
terms of the order A?, eqn. (12) can be solved in terms of a function 

r): f4y=full1+AW,,]. Expanding W, in associated Legendre 
polynomials only one term in the expansion makes a contribution to 


pcos byt Wr, 9, 6)=Uy(r)Po1(cos 9) sin d6+O0(M?4). oe NE) 
In accordance with previous work (Hisenschitz 1949) the function w,, 
has to satisfy the condition 

| OT =e LOL ai Oe) eae ee es (14) 


At small values of r (r<B+3x10-® cm) uy may be replaced by zero 
without a serious error being made. 
Equation (12) is at first solved on the assumption that M=0; in this 


ee. uo =C g(r) +C oh(r) + yrmr®| Q&T 
where C, and C, are constants and the functions g(r) and A(r) are 
solutions of the equation 
r-2(d/dr)r2 du/dr—(d®,/kTdr)(dujdr)—6u/r?=0 . . . . (15) 
g=0 for r=0; h=0 for r+co. It was found that a constant a can be 
Deane so that iim 4 G =a bee tem ee 2 Ses *(L 6) 
=0 


[In eqns. (16), (18), (22), (23) differentiation is indicated by a dash. ] 

In those terms of eqn. (12) which are proportional to M?, a linear 
combination of the functions w,(7) and g(r) is now substituted for the 
unknown w,,(7) so that the solution of (12) which satisfies condition (14) 
is completely specified in terms of known functions 


ray tto IPD) {0x9 + Oph |” GleytaeIHls) Wray ls} (17) 
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where G(r) =(ym/3kT)[ D4’ + (284'/7)—(18®,/77)] 
+ (ug + C399’) (Py! +2B4'/7)+ (Wot C'39)(6P 5/77?) - (48) 
and 


C,=4ymafkT C= ic G(s) fo(s)h(s)s? ds Cg= —| G(s) fo(s)9(s)s* an i 


§5. THe INCREMENT IN VISCOSITY 


The coefficient of viscosity, 7, is equal to 


n=—(N2/24V) | (Oxu/82\yfuWadr . . . . (20) 

and Ay, the increment which is due to the electric field is given by 
Anin=Q ulQe» 15 sees ee 

where 

Qu= | [(etar—ta) x0" + (12MPu9/ 74) — (2M2D,! xo’ to) Th2T?) for? dr 
.. ee ee fess 
Q= | * yf dre 5 ee 

JB 


Expression (21) was evaluated for chloroform. It was assumed that 
the radial distribution f, is similar to that of monatomic liquids having 
a maximum (at r=,), a minimum (at r=R,>R,) and then approaches 
gradually a constant value. A distribution function was plotted so 
that the value of R, was in accordance with the mean distance between 
molecules as derived from the density and that the maximum of fj 
corresponded to 10 nearest neighbours for a molecule. The function ®) 
which is derived from this distribution *s specified in the first row of 
the table. The justification for the use of an almost arbitrary distribution 
is given in the next section. 

No account was taken of the possible effect of the field on the friction 
constant y. Since this constant is closely related to the mean square 
force in thermal equilibrium and the latter is in first approximation 
independent of the field strength it is not likely that, in neglecting the 
field effect on y, an appreciable error was introduced. 


$6. NumericaL Trsts 


The calculated value of 4y/nF? is virtually independent of the exact 
shape of the ®)-curve. This was shown by calculations made with a 
number of different ®-curves characteristics of which are given in the 
table. The viscosities calculated in these tests varied by 1/600 whereas 
the values of 4y/nl? were changed by hardly more than 20%. 

The two approximations to ®) (eqns. (7) and (10)) were compared by 
calculating the function ®,, using the data of the first row of the table. 
Result (fig. 1) shows that the difference is not excessive. The values of 
4n/nF? as obtained with the different approximations have a ratio of 1-33. 
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The calculated relative increment in viscosity is accordingly not very 
sensitive towards the details of the calculation and essentially determined 
by the induced forces of molecular interaction. 


r>R, 
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Ist approximation 2nd approximation 
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§7. COMPARISON WITH EXPERIMENT 


Using the @,-curve as specified in the first row of the table, the values 

M/F=1-10x 10-% e.s.u., -=5-05 and eqn. (10) for Py, it is found that 

An/nF 2=20-'8K 10-8 es... 1 2 ae eeeasy 
comparing with the experimental ‘ viscoelectric constant ° of 
17:0 10-8 e.s.u. for chloroform. Sign and order of magnitude are 
accordingly in agreement with experiment. 

Further experimental results are available for chlorobenzene and 
amyl acetate and do not differ very much from the figure given for 
chloroform. It would appear that the theory is well compatible with 
these measurements although no calculations were made. 

The success of the theory is qualified by the absence of any observable 
effect of the field on the viscosity of benzene, the only non-polar liquid 
investigated by Andrade and Dodd. The corresponding viscoelectric 
constant cannot be larger than 10-® whereas the theory should yield a 
value much smaller than for chloroform but still about 10-8. 

This experimental result suggests an explanation of the field effect on 
viscosity in which the difference of non-polar and polar liquids is 
emphasized and the non-spherical shape of polar molecules plays an 
important part. There can, however, be no gross miscalculation in the 
present estimate of the induced forces and, whatever alternative influences 
may be important, it would be difficult to understand why the induced 
forces fail to show in viscosity. 

Discrepancies would be resolved if there are forces, so far not taken into 
account, which counteract the forces that have been considered. The 
repulsive overlap forces may be of this kind, since the electronic wave 
functions are extended in direction of the field and, correspondingly, 
contracted in perpendicular directions. The exchange integral accounting 
for repulsion should, accordingly, be increased near the induced poles 
and reduced about the equator thus accounting for the behaviour of 
benzene. (The total induced energy and, accordingly, the internal field 
strength are not appreciably affected on account of the short range of 
the repulsive forces.) A corresponding reduction of the viscoelectric 
constant of chloroform would be partly offset by the effect of its 
non-spherical shape so that the above eqn. (24) may be regarded as valid. 
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SUMMARY 


The results of Salzburg, Zwanzig and Kirkwood relating to the molecular 
distribution functions of a one-dimensional fluid are obtained as corollaries 
of a recent contribution by the author and P. M. Mathews (1953) to 
this Magazine. These functions are identified as product densities and 
the problem is considered from the standpoint of stochastic theory. 


$1. INTRODUCTION 


RECENTLY attempts have been made by many authors (see Born and 
Green 1949) to develop a kinetic theory of liquids based on molecular 
distribution functions (they were first referred to as multiform distri- 
bution functions by Born and Green to distinguish them from the 
familiar probability distribution functions). In order to understand the 
nature of the mathematical problems presented by this new theory, 
Salzburg, Zwanzig and Kirkwood (SZK, 1953) studied the case of a 
one-dimensional fluid. In this paper we shall show that a proper inter- 
pretation of these functions from the standpoint of the theory of stochastic 
processes not only rationalizes the method of SZK but helps us to obtain 
their results in the case of a one-dimensional fluid in a very easy manner 
by a new mathematical approach which, it is hoped, will yield fruitful 
results even in more complicated cases. 

In a series of paperst (A, B, C, D) the author formulated and developed 
“the theory of product densities ” in dealing with a stochastic variable 
representing the number of ‘ particles’ or ‘ entities’ distributed in a 
continuous infinity of states. In Paper D a class of product densities 
was considered and we shall now show that the problem of the one-dimen- 
sional fluid if suitably interpreted can be solved as a corollary of the results 
obtained in that paper. For this purpose we shall formulate and solve 
the following problem of random points in a line. 


* Communicated by the Author. roa 

+ Alladi Ramakrishnan, 1950, Proc. Camb. Phil. Soc., 46, 595 ; 1952, Ibid., 

48, 451; 1953, Ibid., 49, 473; Alladi Ramakrishnan and P. M. Mathews, 
Phil. Mag., 44, 1122. 

seen cal hereafter be referred to as Papers A, B, Cand D respectively. 

The notation here corresponds to that used in Paper D, but no attempt is made 

to keep it in conformity with those of A, B and C though we make use of the 


results derived therein. 
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§2. Ranpom Pornts In A LINE 


Consider points randomly distributed on a line represented by the 
parameter ¢ in the following manner. 

The conditional probability that there is no random point between 
i=t, and t=t, (tf, >t,). given that there is a point at t=t,, is ¢(t,—t,) and is 
independent of the distribution of points at parametric values <1. 
It was shown in Paper D that this function uniquely determines the 
statistical distribution of the points on the line. Let z(n,t) be the 
probability that there are n random points between t=0 and ¢ given that 
there is a point at t=0. Using the ‘ method of regeneration points ° the 
author obtained the integral equation 


off 


a(n, t)=— | a(n—1, t—7)h'(r) dr+8(n)o(), «© «© « (1) 
“0 
where ¢'(7)=dd(r7)/dr and 6(n)=1 for n=0 and is zero otherwise. We 
shall now Geitee some corollaries of the above equation : 

If w(n,t) is the probability that the nth point (omitting the point 
at t=0) lies between ¢ and ¢t+d¢ given that at t=0 there is a point. then 
w(n, t) satisfies the integral equation 

t 


w(n, t)= | w(n—1, t—T)w(7r) dr +6(n—1)w(), Sere Fi 


“0 


bo 


) 


where w(7)= w(1, rT) =—¢'(7 
If o(s) and p(n, s) are the eae Transforms* of 4(f) and 7(n, t) then 
it was proved in paper D that 


p(n, 8)=o(8)[1—so(s)|8.> ee ee 
If we define y(n, s) as the L.T. of w(n, t) and write (s) for (1, s), we 
obtain from (2) 
n(n, 8)=[1—so(s)"=[n(s)}". eve Senta) 
From (4) it follows that 


“t ct nt 
a(n =| | oe] eoltsoltety)-. (tty 1) dty dtp sda (6) 
i.e. w(n,t) is the n-fold convolution of the function w(t) where w(t) dt 
is the probability that the ‘first’ point (counted from t=0) occurs 
between ¢ and ¢-+-dt given that a point (uncounted) occurs at t=0. 

We shall call w(t) the ‘ basic function’ of the random distribution of 
points. We shall show that it is w(n, t) and not 7(n, t) that plays a funda- 
mental role in the kinetic theory of a one-dimensional fluid. We note 
that w(n,¢) unlike 7(n, t) is a probability density with respect to ¢ and 


- oO 


| w(n, dt =1 for all n. 
0 


ee ee ee 
* In this paper, L.T, means Laplace Transform. 
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We obtain from (3) and (4) 
sp(n, 8)=y(n, s)—7(n-+1, s) 


: On(n, t 

ie. zit eer Hove iat) sets nee 278) 

This paueeon can also be obtained by elementary probability arguments. 
tee + Gtalk...  t,)* is the product density of degree n of random 


points, given that there is a point at ‘=0, the nature of the random 
distribution we have assumed implies that 


m (bys tay ess b= PO) FP(te—ty) fF Plta—tr—a)> (7) 
provided ¢,,t),...,t, are ordered. As statedin Paper D if t,, t,,... Ags 
are not ordered, f‘/” still exists. f{ by definition is symmetric in the 
variables ¢,,¢,,...,¢,. We can therefore order them and then apply (7). 

By simple reasoning we note that 
fPO= F a(n, t). Pevih et ei aie 
n=1 


Therefore if »s(s) is the L.T. of f(t) 
te ee) 808) 
h(s) ee Th = EO (9) 
For reasons which will be obvious presently we shall call the above 
distribution of points on the line a ‘ free’ distribution with a basic function 
w(t). Stochastic processes associated with such a distribution in the 
special case when w(t)—e~“!A i.e. where 7(n, t) is Poissonian were discussed 
in Paper C. 

Now let us set ourselves the following problem : Given that the N+ 1th 
particle is at =L, what is the nature of the distribution of the N particles 
between 0 and L ? 

For convenience we shall call this conditional distribution ‘ the 
constrained’ distribution with the basic function w(t). For our purpose 
it is sufficient if we are able to determine f\, the product density of 
degree n of the ‘constrained’ distribution. By direct probability 
arguments we obtain 


UNigitias pity 5 L4) 


nr 


5 pol, t,)w(Ms, to—ty) Leg w(Mn; ty —ty—1)O(Mn+ 1> Lia) ; (10) 
aay w(N+1, L) i 
g . 
where »' indicates summation with respect to all m,’s ((=1, 2,...,+1) 


ms 
v n+1 


over all positive integral values of m,(m;40) such that 2’'m,=N-+1. 


t=1 
We shall show in $3 that this is equivalent to the expression for /\ 
obtained by SZK by a tedious process of integration. 
pte Sa ee 
* In this paper f,, denotes a product density of degree n. The superscripts 
are used to distinguish between various product-densities of the same degree. 
For the definition and properties of product densities see Papers A and B. 
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We note two important features characteristic of /'. \ 

(1) fis no longer related in a simple manner to f{? as f “) is related 
to fi. 

(2) Though f@ is uniquely ‘generated’ by w(m,t) the converse is 
not true. In (10), let us replace w(m,t) by another function Q(m, t) 
which ‘ generates ’ the same f‘. We can always write 

Q(m, th=q(m, t)w(m, t). a a ee EL 
Then we find that 
q(m4, ty)q(mo, to—ty) Sipe q(My,; tr—th—1)Q(Mpt 1> L—t,)=q(N +1, L). 


(12) 
n+1 
Since © m,;=N-+1 and L is the sum of the intervals ¢,,t,—t,,..., 
i=l 
L—t,, it follows that 
g(m, t)=exp (um-PAl), an 2 eee 


where » and A are arbitrary constants. We can express (11) in the 
following form : 


Lemma l : 


If Q(m, t) generates f, then exp (um-+At)Q(m, t) also generates f. 
Since w(n, t) is a probability density, if we set 


w(n, t)h=exp [—(un+At)]Q(n, t), «) eee 
( exp [—(un+aAt)]Q(n, t)dt=1, or exp (un)=|- exp (—At)Q(n, t) dt. 
0 0 
(15) 


If Q(n, s) is the L.T. of Q(n, t) from (14) we obtain 
2(n, 8)=exp (un)n(n, 8—A)=[exp (14)n(s—A)]"=[Q(1, 8)"=[.Q(s)]”. 
(16) 
Hence writing 


Q(t) = Q(1, 2b), <b icy al tee nn Oa 
if follows that 


Q(n, =| |. + [ QE )0lte=t)... Olt.) ates aac a) 


“ty—2 


Since w(t) is a probability density 


fo a—| exp [—(u+At)]Q(é) dt=1, i.e. exp (u)= Q(A), sw (19) 
and 
Q2(s+A 
n(3)= on is 22(8)= Q(A)n(s—d). ces ahaa Leo) 
Lemma IT : 


The basic function w(t) can always be obtained from any Q(t) where 
Q(n, t), the n-fold convolute of Q(t) generates the ‘ constrained ’ distribution. 

We note that the basic function contains one arbitrary constant 2 
which we shall call the parameter of the ‘ free ’ distribution. 
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Our fundamental equation (10) for f can be written down directly 
by simple probability arguments only when we deal with the basic 
function. SZK, on the contrary, dealt with a Q(t) which had no pro- 
bability significance and so they were compelled to obtain f@ from 
FY by a process of integration. 

A simple corollary of our eqn. (10) is the following. If w(n, t, L) dt 
is the probability that the nth particle occurs between ¢ and ¢-+dt given 
that the N-+1th lies at L then we obtain 


w(n, thw(N+1—n, LD—t) 


w(n, Up L)= SIN By —, (21) 
and 
(c) 
ee =w(n, t, L)—w(n+1,t,L), . . . (22) 


where 7(n, t, L) is the probability that n particles lie between 0 and ¢ 
given that the N-+1th particle is at LZ. There is of course a particle 
at t=0 and this is not counted. 


§3. APPLICATION TO THE THEORY OF A ONE-DIMENSIONAL FLUID 
SZK assumed an interaction only between the nearest neighbours of a 
one-dimensional fluid. They defined Py(t,, tg,..., ty) dt, dtz...dty as 
the probability that there is a particle between ¢, and t,-+-dt, one between 
t, and t,+dt,,..., and one between ty and ty+dty given that there is 
a particle (uncounted) at {=0 and the N+ 1th particle is at L (this takes 
care of wall effects). They wrote 


Py(ty, ty, .-»» ty) =exp (BA—B[V (ty) + V(t-—ty) 


pe HV (tytn) a VL ty) - + + (23) 
where d 
= L Lh 
exp (—BA)= | |... | exp (PLM) + Mth) 
ee ta, Vl teat diya (24) 


Though it may seem trivial, the author wishes to stress that Py is a 
product density and hence if there are N particles between 0 and L the 
integral with respect to each of the variables ty, ty, ... , ty over the whole 
domain 0 to L should yield NV! and not 1 as it does according to (23). 
Py should be identical with f\? which according to (10) should be 


given by 


Py=fPlty tas ---> ty) 
__ exp (—B[ Vt) +VG.—th)+- . EV (tw—ty- 1) + Absa (25) 
Q(N +1, L) 
where 
gin+1.B)=[ fof exp (et )+ Mt) 
‘ Drs) tn—1 


+. . ot V(ty—ty—1) + V(L—ty)]) dt, dts ears dt x. . . (26) 
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In spite of this, SZK obtained correct expressions for f© (which 
they denoted by p) because they consistently integrated expressions. 
involving the numerator of their Py over the whole domain thus intro- 
ducing a factor NV! in the numerator, which annulled the defect in their 
definition of Py. We will rationalize their process in §4. Comparing 
the corrected expression for P= f? with (10) we find 


QO) = exp[ SVQ) =] = ee 

Thus the basic function is given by 
o(t) exp (—[us-A-E BVI) ev 
sates of (=| exp (-TM-LBVD die ee ee 


We note that « is a function of A and hence our expression for the basic 
function w contains one arbitrary constant which is consistent with the 
statement that the potential V contains an arbitrary constant. Thus 
the expression for f\ given by SZK is identical with (10). The 
corrected expression is Py = f= p™.* 
Now let us consider the following potential V : 
V(é)==00, 0<i<a 


V(t)=0, t>a «ge Feed oy ee 
Le; exp [— V (t)]=A(t—a), : (31) 
where H(x) is the Heaviside unit function, i.e. H(x~)=1 if «<>0, H(x)=0 
if x<0. 
aS w(t) =A (t—a)A exp [—A(t—a)], —p=logA+aA. . (32) 
— na n—-1 
w(n, t)=H(t—na) exp [—A(t—na) | a Piette ee CH 
is given by applying (10). In particular 
fl N! (t—na)"-1 (L—t—N+1—na)*—- 
Ni) —foOnh— > A RINGS eR ee eee 
PO=lrO= 2 Gop! W—n)! (LN ay 
x H(t—na)H(L—t—N-+l—na). . . . . . (84) 


§4. InreGRaTION Process or SZK 
We shall now discuss the integration process of SZK. For this purpose 
we have to recall another species of distribution function—the Janossy 
density defined in Paper B. Adapting the definition of Janossy density 
to our problem we say J,,(t;,t,...,¢,) dt, dt,...dt, is the probability 
that there is one particle in dt,, one in dt, ... , and one in dt, and none 
elsewhere. We quote here a result proved in paper B connecting product 
densities and the Janossy densities : 
oe) il DoE yf, L 

Frlta by, = Tenet ee (m—n)! |. if Set | J in(ts, ty, Sw ca em) 

UO i each | RR ee le EN 


* If we use the SZK definition of P,,, p™ will not be identical with Py, and 
this is illogical. ; N 


— 
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If we are given that the N+1th particle occurs at L the Janossy densities 
of all orders except m=N vanish and 


' 1 L pL -L 
(n) — fo) 2 
Pp wv = bs (i, to, er) b= (N—n)! i fe erie. id | Seta. ty, a) ty) 


Sit amrkigy ie aia cwittine aaa ete (36) 


We also note that JQ=/® and this is identical with our corrected 
definition for Py. The most interesting feature of J® is that it is 
symmetric in the N variables ¢,, ty, ... , ty and exists even when they are 
unordered but can be explicitly written down according to (25) only 
after we order them. Hence the computation of the multiple integral is 
complicated and is done in the following manner. 

We first order the N—n variables t,,,, t,.9,.--,¢y not concerning 
ourselves with their order relative to the other n variables Beg UP Resta chat 
We rewrite eqn. (36) using the symmetry property of f@ as 
FO (Ea tay +» tn) 

rL ph 
= | 

Let us now assume ¢,,t,,...,t, are ordered. J cannot yet be 
explicitly written down unless the order of ¢,, ¢,,...,#, relative to the 
other N—n ordered variables ¢,,.1,t,19,...,ty is known. The relative 
order changes in the process of integration and so to evaluate the integral 
(37), SZK adopted the ingenious device of writing (37) in the form 


oe elle FWY, L 
OE Ree a | [ ne | pai feel eaves oy thee Ty): 


ne 


L 
he | Ohta Ua wae eee dip les ete eee (37) 
tt ~tny—1 


2 mae 4 aay 4 
OTe 13)0(T 4, — be) Olt tn) O71 OTe Oty. - = » (38) 
+ denotes summation over all values of k,,k,,...,k, subject to the 
k 
condition k,<k,<...<k,. Since according to (25), JY is given in 


terms of Q(r)=exp [—fV(z)], this process yields our fundamental eqn. (10) 
but with Q(n, t) in the place of w(n, t). Though SZK used the uncorrected 
expression for J?= Py in (44) they were led to the correct result 
since in (44) they integrated with respect to each of the V variables over 
the whole domain 0 to L thus introducing the factor N ! which annulled 
the defect in their definition of Py. 


§5. Asymproric Forms or f\? = p\™ 
If N>1; Lt, a (34) takes the asymptotic form* 


ve 
0 ox ae aN = pall oN 
pY(t)h~ 2 H(t—na) ied ee? [S| ; (39) 


to (n—1)! L—Na 1 ING 


* We shall not bother ourselves with a formal derivation of the asymptotic 
form which can be obtained in just the same manner as a Poisson approximation 
from a binomial distribution. 
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We can also write the asymptotic expression for p as 
pM (ty, to, OPC) tn) es p')(ty)p\ (t2—th) oes peal sited + 59.) : u (40) 
provided N>1; L>t,, to—t, ts—ty,.--,ty—ty-1 and L—t,, >a. 
Comparing (39) and (40) with the corresponding expressions in the case 
of a ‘ free’ distribution we are led to the following result : 


Lemma IT] : 

In the case of a potential of width a and infinite height if N>1 and L>a, 
as long as we are dealing with the distribution of particles in a domain 0 to t 
where L>t, we can replace the ‘ constrained’ distribution by a * free’ 
distribution with the parameter A4=p=N/(L—Na). 

We therefore write, using (22) and (33) 
tr-1 


rt—na 
P An exp (—At) Ca 


(eel 


n(n, t, L)= | 


N 
L—Na’ 

SZK obtained the asymptotic expression for p'?)(t,, t)) in the case when 
t,>a, but t,—t, is comparable to a. They wrote p(t,) ~ N/L for 
t,, L—t,>a,t,—t,. The author is not in agreement with SZK in regard 
to this value of p‘)(t,) and wishes to point out that p')(t,) ~ N/L—Na 
for t,, L—t,>a. A comparison between the asymptotic forms of SZK 
and that of the author will be made in a later paper in which he also 
proposes to discuss the molecular distribution in the case of a square 
well potential defined by 


A"+1 exp (—At) — dt, A= (41) 


“80 


V(t)= 00, 0<t<a 
Vij==e,. axt<b 
Viij=0; t>b Ltd 1 Se ee) 


by using the method of basic functions. In that paper it will be shown 
that if the basic function w(t) is discontinuous at a finite number of 
points, it leads to difference differential equations which can be com- 
pletely solved. The problem of a one-dimensional gas with potential 
functions defined by (42) was proposed and discussed recently be Sells, 
Harris and Guth (1953). 


$6. Some Remarks ON PRopucT DEnsITIES 


Since molecular distribution functions have been identified as product 
densities, it is relevant to note the following properties of product density 
functions in addition to those stated and proved in Papers A to D. 

(1) Product densities of any degree are defined with respect to a 
continuum, which in general can be n-dimensional. If a product density 
is integrated with respect to 7 dimensions it is still a product-density with 
respect to the continuum of n—r dimensions. 


— 
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(2) Taking in general a product density of degree Nias /j(t,, fas». 5 t,)* 
Jet us ask the question : What is the effect of integrating f,, with respect 
to only one variable? Writing 


if VAUD ty, ee eg t,) dt, dt. ok ae dt. =F, ita» to, ee ey ty 1) dt, dts wave dt, 1 
aa (3 


(43) 
this represents the expectation value of the product of the number of 
particles in the range of integration of ¢, and the number of particles in 
dt, dts,...,dt,_,. Then R,_, is no longer a product-—density and 
RK, Ut, dt, . . . dt,_, no longer a probability magnitude. 

In particuiar when we have a total number N and only N particles 
in the entire range, 

RY? (ty, bg, Me eats t,-J=(N—A+1) M(t; bas siises9 ty-1)> (44) 
if f;, is integrated over the whole range of t,, since the number of particles 
outside the infinitesimal ranges dt,, dt,,..., dt,_, is N—h-+-1 conditional 
upon the existence of particles in dt,, dt,,...,dt,_, and so it is not a 
random variable. From (44) we immediately get 


1 
(Ents, icae, t =a—m/! | soe SP ie bas eens tie) 
Pits 2 n) (N—h)! pre es oe 912, N 


eel Olt as neni, ee ee cae 20) 
Born and Green state that (44) is obvious but it is not clear why it is 
so unless we interpret their functions as product densities. Since such 
functions were not so identified and distinguished from probability 
densities, in the case of the ‘ distribution functions of many electron 
systems ’, Mott and Sneddon (1950) had to appeal to an artificial rule of 
integration. They stated that to obtain f\ from f{ we should 
integrate fQ with respect to t,+1,t,+2,..-,¢y and the “ integration 
must be carried out in such a way that each configuration must appear 
only once ; that is, dt,+ 1 dt,+....dty must be replaced by 
dirt Uap we OLN ss 
(N—h)! 
This argument is not correct and must be replaced by the mathematical 
arguments outlined above. 
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Produced in Steel by 340 Mev Proton Bombardment 
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Londonderry Laboratory for Radiochemistry, University of Durham 
and 


K. I. MAYNE 
Clarendon Laboratory, Oxford* 


[Received January 8, 1954] 


SUMMARY 


The small amounts of 4He, ?He, and #H, produced in a block of steel 
bombarded with 340 Mev protons have been isolated. Measurements. 
have been made of the total helium by a gas-volumetric method, of the 
tritium by counting of the associated f-particles, and of the *He : *He 
ratio by mass-spectrometry. The ratio *He : *He : 3H has been found to 
be 1-00 : 0-048 : 0-078 ; this result is discussed in the light of Le Couteur’s 
theory of nuclear evaporation. 


INVESTIGATIONS of the nuclear evaporation processes induced by 
bombardment with particles of very high energy have, so far, mostly 
been made by one or other of two complementary techniques: by 
radiochemical studies of the yields of radioactive spallation products 
remaining after the completion of the evaporation process, and by the 
examination of the tracks produced in photographic plates by the charged 
particles which have been emitted. 

The correlation of the resulting data is rendered difficult by two facts. 
In the first place, ‘evaporation stars’ in photographic plates must, in 
general, have originated in one of the elements (H, C, N, O, Br, Ag) 
naturally present in the emulsions. This immediately restricts the 
range of possible comparisons and, further, makes it difficult to obtain 
an unbiased selection of ‘stars’ without introducing some uncertainty 
concerning the origin of the events. It therefore seemed desirable to 
make a comparison of the light particles evaporated with the spallation 
products from the same material under similar conditions of bombardment. 
Secondly, although it is comparatively easy to distinguish between 
singly and doubly charged particles, mass discrimination is by no means 
so readily achieved. Most of the attempts in this direction have involved 


* Communicated by the Authors. 
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some inherent bias, introduced, for example, by the necessity to select: 
tracks of length sufficient for the measurement of the mean deviation 
from a rectilinear path (see, for example, Holtebekk, Isachsen and 
Sdrensen 1953). 

In the present work we have made a direct measurement of the total 
yields of three of the evaporation particles (*He, 3He and 3H) produced. 
in a sample of steel which had been bombarded with 13 yA-hr. of 340 Mev 
protons in the Berkeley synchro-cyclotron: 

Total helium has been measured gas-volumetrically in an apparatus 
similar to that recently described by Gliickauf (Gliickauf 1946, Chackett, 
Reasbeck and Wilson 1953), tritium has been estimated by gas counting 
of the f-particles emitted (Allen and Ruben 1942, Libby and Barter 
1942), and the helium isotope ratio determined by mass-spectrometry 
in a high resolution Nier-type instrument (Mayne 1954). Helium was 
extracted by solution of the steel in aqueous copper potassium chloride 
solution ; this resulted in an excessive dilution of the tritium with 
ordinary hydrogen, and the release of the tritium was therefore achieved 
in the dry way. Vacuum fusion, and ignition of the steel in an atmosphere 
of chlorine were both used, and gave substantially the same result. 
Since it was thus not feasible to measure helium and tritium on the 
same sample, and since the irradiation of the block was very non-uniform, 
the results have been inter-related by measurement, in a counting system 
of reproducible geometry, of the gross radioactivity of the spallation 
products in the sample. 

From that part of the steel on which the cyclotron beam had been 
directly incident we find (from pairs of concordant measurements) 
4He : He : 3H: : 1-00: 0-048: 0-078; from another part, which was 
probably irradiated mainly by scattered particles (of rather lower energy) 
somewhat smaller ratios were found: *He:?He:?H:: 1-00: 0-028: 0-038. 
Correction has been made in each case for decay of 3H in the interval 
between irradiation and measurement. 

The correlation of these results with the predictions of nuclear evaporation 
theory (Le Couteur 1950, 1952) is rendered difficult by the fact that the 
theoretical results are necessarily expressed in terms of the energy of 
excitation of the nuclei. Protons of 340 Mev are sufficiently energetic 
to excite meson production, and ‘knocked-on’ nucleons may also be 
emitted before the true evaporation process commences; it is thus 
likely that in our sample the average excitation energy was considerably 
less than the full 340 Mev. 

We have attempted to make allowance for this uncertainty by utilizing 
the data of Rudstam, Stevenson and Folger (1952) on the yields of 
spallation products from iron bombarded with 340 Mev protons. From 
their yields (with appropriate interpolation where necessary) one can 
deduce cross sections for events leading to the evaporation of various. 
numbers of nucleons. Le Couteur (1950) has given graphical representations 
of the average yields of various evaporation products, as functions of the 


Z2H2 
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energy of excitation of the nucleus. From his curves, one can deduce 
relationships between the total number of nucleons lost by evaporation 
in a given event, and the average numbers of particles of various sorts 
constituting this total. Combining this with the cross section data of 
Rudstam et al., one can then arrive at overall weighted average yields 
for the types of particle with which we are concerned. In this way we 
have deduced ‘He ; #He : 3H: : 1-00: 0-068: 0-178. In view of the 
assumptions involved, the agreement with our experimental results 
is not unsatisfactory. Our experimental value for #H :4He is rather 
low compared to the calculated value; this may arise from the fact 
that Le Couteur’s calculations (which refer specifically to a hypothetical 
atom of mass 100, excited so as to produce at least seven charged 
evaporation particles) have here been extrapolated to cover an iron 
nucleus (4=56) at moderate levels of excitation. It is to be expected 
that the change from A100 to A=56, with the consequent reduction 
in potential barrier, would tend to reduce the yield of singly charged 
particles relative to doubly charged fragments. 


Our thanks are due to Professor F. A. Paneth, F.R.S., for his 
interest, to Professor 8. F. Singer (University of Maryland) who suggested 
the experiment, to Professor W. H. Barkas (University of California 
Radiation Laboratory) who kindly undertook the bombardment of the 
steel in the synchro-cyclotron, and to Messrs. Imperial Chemical 
Industries, Ltd. who have generously made available to us some of the 
equipment used in this investigation. 

One of us (K. I. M.) is indebted to the Pressed Steel Company, Ltd., 
for financial support ; another (G. W.) wishes to acknowledge the receipt 
of an award from the Ministry of Education, during the tenure of which 
this work was carried out. 
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Emulsions* 
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SUMMARY 


In an examination of two-particle events recorded in photographic 
emulsions, four examples have been observed of a type similar to the 
decay of A°-particles in Wilson chambers. In three of them, the measured 
energy-release in the decay is ~40 Mev; in the fourth example it is 
~70 Mev. An event has also been observed, which appears to correspond 
to the decay in flight of a 6°-particle, emitted from a nuclear disintegration 
of type 10--7p. One of the secondary particles arising from the decay is a 

“negative 7-meson. The measured Q value is 202+11 mey. 


§1. Tor ENERGY RELEASE IN THE Decay oF A°-PaRTICLES 


It is now well known, from experiments with Wilson chambers, that a 
two body decay exists of the type :— 
Ao > pta-+Q. 

The available evidence has been summarized by Rochester and Butler 
(1953), and some of the more recent measurements are contained in the 
Bagnéres Report (1953). Values of @ ranging from 34 Mev to 45 Mev have 
been reported by various investigators. The Pasadena group (Leighton 
et al. 1953) suggest that a second Q value of ~75 Mev also exists. We shall 
here describe some determinations of the energy release Q in decays of the 
above type observed in photographic emulsions. 

In this laboratory, observers have been recording two-prong events, 
collectively referred to as 2+0 stars. These events belong to one of the 
following categories : 

(a) elastic or inelastic scattering suffered by a charged particle ; . 

(b) neutron induced nuclear disintegrations resulting in the emission of 

two charged particles ; 

(c) one-prong stars produced by a charged particle ; and 

(d) decays of various types of V-particles. 


* The nomenclature used in this paper is that suggested by Amaldi et al. 


1954). 
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In order to isolate events of type (d) from amongst the various 2--0 
events recorded, the following two criteria were used : 

(a) Only those 2-0 events were selected in which it was possible to 
show, by combined measurements of grain-density and multiple scattering, 
that one of the tracks was due to an outgoing L-meson ; 

(6) no events were accepted in which there was a recoil track or spur 
of length >0-5 p, or an associated electron or 5-ray at the apex of the fork. 
Amongst the large number of 2-prong events examined, 4 were found 
which satisfied the above criteria. Details of these events are given in 
table 1. 

Table 1 


Track 
length in g* pp Identity ie E @ Q 
microns 
1 CM (a) 800 2-23+0-074 60+19 L-meson 102+19 33+11 124° 53+13 
2 (b) 1900 3-28--0-:068 101417 Proton 321+30 54+10 


Event Expo- 
sure 


) 
> OM (a) 1000 2-73+0-066 54+12 L-meson 96+13 30+ 7 
7 : (b) 1400 3-88+0-083 95+19 Proton 304432 49+10 127° 43+11 
3 eT A 24 (a) 1800 3:15+0-063 25+ 6 L-meson 62+ 7 13+ 3 
“—“* (6b) 300 3-8 +0-128 Proton 137 oo 
4 OS (a) 4880 3-2 +0-07 m--meson 68+ 1 16+ 0-5 
(b) 2800 1:98--0-05 433-+87 Proton 712+88 242+ 53 135° 80+10 


CM: = Ilford G5 emulsions, 600 thick, exposed for 50 h, between 30000 and 42 000 ft., 
with 14 cm of iron on the top and 9 em of iron on all sides, in a Comet jet airliner 
flying between London and Singapore. 

HA3¢4: Ilford G5 emulsions, 400 » thick, exposed for 5h at ~65 000 ft., under 10 em of lead. 

ACIS) A block of Ilford G5 stripped emulsions, each 600 u thick, exposed for 25h between 
30 000 and 42 000 ft. in a Comet jet airliner flying between London and Singapore. 


Values of g* given in col. 4 are the normalized grain-densities. Tracks 
of high energy electron pairs were used for the normalization and for 
determining the variation of development with depth in the emulsion. In 
event 4, track (a) is due to a negative 7-meson, which after traversing 
4-88 mm in the emulsion was brought to rest and produced a 2-prong star. 
In the other examples, it cannot be stated whether the L-meson is a z- 
or a -meson ; nor can the sign of the charge of the meson be determined. 
The identities given in col. 6 are based on the grain-density and scattering 
measurements. In event 3, no scattering measurements could be 
carried out on track (b). It has been attributed to a proton, however, 
since it could be seen visually, that for a measured value of g*—3-8, the 
track was not sufficiently scattered for it to have been produced by an 
L-meson. The space angle ® between the two tracks in each of the 
examples is given in col. 9. In all the events the measurements favour 
the assumption that the L-meson is travelling away from the event, but 
in some the evidence is not decisive. In the case of event 3, however. it 
can be stated unambiguously that track (a) was produced by an outgoing 
m-meson. The @ values given in col. 10 have been calculated on the 
assumption that the L-meson is a 7-meson. 

The values of the momentum, p, and energy, H, given in cols. 6 and 7 
of table 1 are based on scattering measurements. It is obvious, however, 


Decay of Neutral V-Particles in Photographic Emulsions 415 


that such measurements on limited lengths of track, are subject to large 
Statistical fluctuations. In order to determine the @ values more accur- 
ately, an alternative procedure was therefore employed. In each case, it 
was assumed that track (a) was due to a7-meson, and track (b) to a proton. 
Then, using the values of g* given in col. 3 of table 1, values of pB were 
deduced from the best g* vs. p8 curves which had been obtained experi- 
mentally for protons and 7-mesons in these emulsions. The values of 
p, E and Q, derived from these values of pp are shown in table 2. Photo- 
micrographs of events 1, 2 and 3 are shown in Plate 10, 


Table 2 
Event Particle p E ® Q* 
(a) 7-meson. 88+ 5 DOs B 
I roten 338-411 Hen ce ee oe 
> (a) 7-meson 80 + 4 DO ae 2 127° 4] +2 
~ (b) Proton 260+ 5 352 2 r, riba 
(a) 7-meson 68+ 3 16+ 2 Qr70 Qn 
: (b) Proton 285 +28 47+ 8 Se AE as 
(a) 7-meson 68+ 1 16+ 0:5 50 70) 
4 6) Proton 687420 226412 °° ue 


*The errors shown are standard deviations; they are statistical errors 
associated with the grain density and angle measurements. 


Since the direction of motion of track (6) in each event could not be 
determined unambiguously, we have also considered the following 
possibility. The events were assumed to represent the decay in flight of 
charged hyperons, the decay involving the emission of a charged L-meson, 
and a neutron. The velocity of the hyperon at the point of decay was 
-obtained from the grain density of track (b) of each event. The Q values 
obtained from the four events for a decay of the above type, were: (1) 
23 mev, (2) 17 Mev, (3) 11 Mev and (4) 23 Mev. These values are incon- 
sistent with other determinations of the energy release in the decay of 
charged hyperons. 

A summary of the data presented in this section was given by the 
author at the International Conference on the Cosmic Radiation held at 
Bagnéres-de-Bigorre (July 1953). At the Conference, a fifth event had 
also been reported, in which both the secondaries of the V-event came to 
rest in the emulsion and were identified as a proton and a 7 ~-meson 
respectively ; the Q value obtained in this example was ~10 Mev. 
‘A closer examination of the event suggests, however, that it represents a 
nuclear interaction and not the decay of a A°-particle. In view of this, 
we have excluded it from our statistics. 


§2. Drcay oF A 0°-PARTICLE 
A star of type 10+ 7p was found recently in an emulsion block consisting 
of 46, 600, thick, Ilford G5 stripped emulsions exposed at 75 000 ft. for 
~4hoursin England. Ata distance of ~100u from this disintegration, a 
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two-particle event was found (see Plate 11), which appears to represent 
the decay in flight of a 6°-particle. The details of this V-event are as. 
follows : 

(a) One of the emerging particles was a negative 7-meson, for it came to 
the end of its range in the emulsion after travelling a distance of 3-1 mm 
and produced a one-prong o-star. 


(b) The other track has a total length of 54mm in the stack and has. 
been observed in eleven emulsions ; the minimum length in any single 
emulsion was 4mm. The mean blob-density on the track was determined. 
to be b¥=0-928-+0-019 and multiple scattering measurements yielded a 
value of pB=412-21 mev/c. Individual measurements, carried out on 
each of the sections of the track in the eleven emulsions, were found to be 
consistent with one another. From these values it may be seen that the 
track must have been produced by an L-meson. It is not possible to. 
determine whether the L-meson involved is a z- or a y-particle, nor can its. 
direction of motion be established. 


The simplest explanation would be that the V-event represents the 
inelastic scattering of a negative 7-meson. This possibility appears to be 
very unlikely, however, on the following grounds : 


(a) There is no sign of any recoil or electron track at the apex of the 
fork. If the event were due to an inelastic scattering, the 7-meson has. 
been degraded in energy from 340 Mev to 12 Mev inthe process. That such 
an energy loss should take place without any visible trace, appears to be 
highly improbable. 

(6) The two tracks of the V-event and the centre of the disintegration 
10+-7p are co-planar to within ~2°. 

In view of these features, we suggest that the event is an example of the 
decay in flight of a 6°-particle, the latter probably having been emitted 
from the star 10+-7p. 

The event has been analysed on the basis of this assumption, and the 
relevant details of the measurements are given below. 


Angle between secondary mesons : = 146°, 39'+30° 
Angle ®_ between 7~-meson and assumed direction of 6° =142°, 2'’+30' 
Angle ®, between fast light meson and assumed direction 

of 99 = 4°, 54’-+15" 


(It is assumed that the direction of the @°-particle is given by the line 
joining the apex of the V-event and the centre of the disintegration 
10-+- 7p.) 


Range of negative 7-meson =3128y 
Energy of negative 7-meson =12-2+0-2 mev 
Momentum of negative 7-meson =59-6+0-41 Mev/c 


Transverse momentum of negative 7-meson =37 -+0-5 Mev/c 
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If the fast particle If the fast particle 


is a 7-meson is a p«-meson 
Energy 342 +26 Mev 358 +27 Mev 
Momentum 461427 mev/c 452 +28 mev/c 
Transverse momentum 39+ 4 Mev/c 38:5+ 4 Mev/c 
@ value 202+ 6 Mev 220 + 6 Mev 
Mass of 0° 950 m, 921 m, 


In calculating the energy of the fast particle, account has been taken of the 
loss of energy by ionization along the track. The error quoted for the 
energy of the negative 7-meson is solely that due to straggling. 

The energy release Q involved in the decay 

0° +a7t++n-+Q 

is thus determined to be 202+11 Mev. This value may be compared with 
that of 214 mev quoted by Thompson et al. (1953). 

The error of +6mev has been increased to +11 Mev to allow for 
uncertainties (~5°%,) in the value of the scattering constant employed for 
the estimation of pS, and uncertainties in the thickness of the emulsions. 
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SUMMARY 


Two examples of charged hyperons have been observed in photographic 
emulsions. One decays in flight, and its transformation may be re- 
presented Y* + L*+N°+Q, the release of energy, Q being consistent 
with values obtained in other laboratories. In the other example, the 
secondary is a fast proton. This event may represent either the decay 
of a new type of hyperon or the nuclear capture of a negative hyperon. 


§ 1. INTRODUCTION 


At the Royal Society Conference in January 1953, the first decisive 
evidence from experiments with photographic emulsions for the existence 
of a charged hyperon was presented (Bonetti e¢ al. 1953 a). At the same 
time, and by independent observations with Wilson expansion chambers, 
York e¢ al. (1953) obtained similar evidence. These workers observed 
the decay in flight of an unstable particle produced in a nuclear disinte- 
gration, where the secondary particle had a mass not significantly different 
from that of the proton. Several examples observed in photographic 
emulsions and expansion chambers have since been reported from other 
laboratories. In some cases the hyperon is observed to decay in flight 
whilst in others it reaches the end of its range. The limited evidence 
hitherto available has been consistent with the view that we are dealing 
with the decay of particles of a single type, particles of which the 
transformations may be represented 


Y+ +n+7*+4+Q (a) 
and. Yt > p+79+Q. (b) 
For those examples observed in photographic emulsions, only one can 
be interpreted in terms of (6); all the others may be represented by (a). 
The Q-values for the two modes of decay appear to be equal within the 


present wide limits of error, Q ~ 110-130 Mev, and the assumed two-body 
decay rests upon this evidence (Bonetti et al. 1953 b) 


eee 

* The nomenclature used in this paper follows the system set out by Amaldi 
et al. (1954). 

} Communicated by Professor C. F. Powell, F.R.S. 

t Note added in proof.—A recent event has been observed by Danysz (private 
communication), in which a charged hyperon, emitted from a nuclear disinte- 
gration, decays in flight, the charged secondary particle being identified as a 
proton. The calculated @-value for the decay, according to (6), was 
Q=12018) mev. . 
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The cascade decays observed in expansion chambers (Armenteros 
et al. 1952, Anderson et al. 1953) may also represent the decay in flight 
of a charged hyperon, into a light meson and a neutral particle, either 
Ae or 0°, but this interpretation remains to be confirmed. In each case, 
it depends on an observed coplanarity of one V event of a pair and the 
apex of the other. Whilst strongly suggestive, these observations do not 
prove that the neutral particle results from the decay of the charged 
particle. They might both result from a nuclear interaction in the matter 
above the chamber—a form of pair production. (See for example Fowler 
et al. 1953 b.) 

In all examples of this type, reported hitherto, only lower limits can 
be placed on the mass of the primary particle. Of the primary particles 
of such events, observed in magnetic fields, all are negatively charged. 

Two examples of slow charged hyperons (Y(B)1, Y(B)2) have recently 
been observed in stripped emulsions in this laboratory. In event Y(B)I1, 
shown in Plate 12, a charged hyperon appears to decay in flignt. The 
secondary, although steeply dipping, can be traced through eleven 
emulsions and can be identified as a light meson. In the other event, 
shown in Plate 13, the hyperon comes to rest. The secondary particle 
in this event is also steeply dipping, but can be traced through sixteen 
emulsions and can be identified as a proton. 

These events are described in §§ 2 and 3 and are discussed in relation 
to other charged hyperons in § 4. 


§2. Event Y(B)1 

A photomicrograph of this event is shown in Plate 12. The hyperon 
is ejected from a star of type 26+ 14p and travels 4-72 mm in one emulsion 
before decaying in flight. All other tracks from the star have been 
traced either to rest, or until they left the emulsion stack. All the 
26 heavily ionizing and 10 of the shower particles have been identified 
as protons, deuterons, «-particles or 7-mesons. The remaining 4 fast 
particles could not be identified. The residual range of the hyperon 
at the instant of decay, its mean velocity during flight, and its mass 
have been determined by observations on the variation of gap-density 
and scattering along the track. These results are summarized in table 1. 


Table 1. Measurements on the Primary of Y(B)1 


Mass of Primary Residual range at 


Method tn, decay B=v]c 
Scattering : 2810 +600 6-3+1-:0 mm 0-254 +0-010 
Gap-density : 2670 +600 6-5 +1-0 mm 0-265 +0-010 


As a basis for further calculation, a mass value of 2750 m, and a velocity 
B=0-26 have been assumed. The results thus obtained are not very 
sensitive to small changes in either of these quantities. The proper 
time of flight of the primary was thus found to be 5-8 x 10~™ sec. 

The secondary is directed backwards, in the laboratory system, at a 
true angle of 34:5° to the direction of motion of the primary. Because 
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the track of the secondary is very steeply inclined to the plane of the 
emulsion, an accurate estimate of the energy cannot be obtained from 
scattering measurements. The track can be traced through eleven 
emulsion strips until it leaves the stack, and the mean angle of scattering 
over the entire length is <=0-34°/100u. Blob counts are more reliable.. 
and the mean normalized blob density in the eleven strips is b*= 1-34 +-0-02. 
This value and the observed «, indicates that the secondary particle was. 
almost certainly a light meson. From the value of the blob density, 
and the relation between b* and pf already established for these plates, 
values of p8, assuming the secondary to be either a 7- or u-meson, can be 
deduced. These values are pB=87-+3Mev/c for a m-meson; and 
pB—=66+2 mey/c for a p-meson. The observed energy of the meson. 
allows us to exlude the possibility that the event represents a nuclear 
interaction involving, for example, a slow deuteron, for the parent 
particle has not sufficient kinetic energy. 

Assuming that the event represents the decay in flight of a charged 
hyperon, and that the transformation is into two particles only, the decay 
may be represented 

Xe Lee, 

where L+ denotes a light meson, and N° a neutral particle. @-values 
and mass values for the hyperon, have been calculated for all possible 
combinations, assuming L to be a z- or u-meson, and N® to be one of the 
following neutral particles already known: A® (mass 2180), n (mass. 
1838), B° (mass ~ 1400, Fowler et al. 1953 a), 0° (mass 970) and 7° (mass 
263). These values of Q and My are set out in table 2. If the hyperon 
is assumed to be an excited nucleon, and if it is accepted that nucleons 
are conserved in these processes, then the possibilities that N°—B®, 6° 
or 7° may be immediately excluded. However, if the hyperon is a 
created particle, then all possibilities must be allowed. 


Table 2 
(a) @-values (Mev) 

Secondary \ Neutral Particle 51h n B® 6° a 
7 ; 98 100 106 115 172 
P : 72 73 76 81 119 

(6) Mass Values of Primary (m,) 

Secondary \ Neutral Particle 43 n Be 6° 77° 
7 2645 2306 1883 1468 73 
m 2528 2187 1756 1336 703 


From these figures and the direct mass-measurements on the primary, 
it seems very unlikely that the neutral secondary is less massive than the 
proton. The most probable Q-value is thus ~ 100 Mev. It is, however. 
not possible, with the present information, to decide whether the neutral 
particle is a neutron or a neutral hyperon. If the latter possibility is 
the correct one then the event could represent a process similar to the 
‘cascade decays ’. 
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§3. Event Y(B)2 


In this event (Pl. 13), the hyperon was ejected from a star of type 
31-++3n and came to rest after a path of 4:64 mm in three emulsions. 
All other tracks except the shower particles, have been followed from the 
star and identified as protons, deuterons, «-particles, or 7-mesons. From 
the end of the hyperon track, a second charged particle emerges, with a 
track steeply inclined to the plane of the emulsion, which leaves the stack 
after 13-2mm, having passed through sixteen emulsions. The measure- 
ments on the primary are summarized in table 3. 


Table 3. Mass Measurements on Primary of Y(B)2 


Method : Range-scattering Grain-density-range 
Mass : 2680 +350 m, 2200 +350 m- 
Fig. 


TOTAL NUMBER OF GAPS 


20 


RANGE (M.M) 
Total number of gaps vs. range, from the end of the track. Comparison curves 
are given for deuterons, protons and z- and p-mesons; the curve for the 
hyperon is that denoted by J2. 


In making the range-scattering measurements, three different sets of 
cell sizes were employed in conjunction with the constant sagitta method 
(Biswas et al. 1953). The value quoted is the mean of the three deter- 
minations. Gap measurements all indicate a mass between that of the 
proton and the deuteron, but the experimental conditions render the 
observations unreliable and it is therefore not possible to give a definite 
mass value. A plot of the total number of gaps as a function of range, 
is shown in fig. I, which includes curves for z- and 7-meson, protons and 
deuterons, for comparison. 

As in event Y(B)1, the track of the secondary particle is very steep. 
Tn the first three plates counting from that in which the hyperon comes 
+o rest, the secondary has a mean grain-density of 29-6 grains/50 pu. 
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Scattering measurements give a value of a—0-25°/100 4. These measure- 
ments show that the secondary particle was considerably heavier than 
the 7-meson; and taking account of the fact that the value of « is an 
upper limit, the track being steeply dipping, they are consistent with the 
assumption that it was a proton. If the secondary particle were not a 
proton, but a lighter particle, then in the length observed in the emulsions. 
(13-2mm) a significant change of grain density should be observed. 
No such change has been found. The grain density and gap measure- 
ments have therefore been employed to estimate the total range of the 
proton, and a value R=3-0-+0-8 cm has thus been obtained. This. 
corresponds to a proton energy of 100+ 15 Mey. 

We may consider two possible explanations for this event. If it 
represents the decay of a charged hyperon 


¥+ > pt+m+Q, 


then it is easily shown that the energy release, Q, is of the order of 
430 Mev, a value considerably higher than any previously reported for 


any similar event interpreted as due to the decay of a charged hyperon. 
An alternative explanation is that the event represents the nuclear 
capture of a negative hyperon, the proton being the only visible product 
of the disintegration. 
§ 4. Discussion 


A recent paper by Bonetti et al. (1953 b) summarizes the results. 
available until the beginning of November. Since then two hyperons 
have been observed by King ef al. (1953) ; one of these hyperons comes to 
rest and has a fast 7-meson secondary. Considering two possibilities for 
the neutral secondary for a two body decay, viz. A° and a neutron, they 
obtain Q values of 101 and 103 Mev. Event Y(B)1 is in good agreement. 
with this. The events listed by Bonetti e¢ al. have slightly higher Q 
values, but within the accuracy of all the measurements, limited by 
steepness or short track length of the secondaries, all of the decays so far 
observed are not mutually inconsistent. In appearance event Y(B)1 
is very similar to that observed in Padua (Ceccarelli and Merlin 1953) 
which event has a slightly larger Q value. 

Observations on cascade decays may be taken as supporting the 
identification of the neutral secondary, in the charged hyperon decay, 
asa /°, but the weight of these results is not yet very large.* In this 


* Note added in proof.—At the meeting of the American Physical Society 
in New York (January, 1954), Cowan presented the results of a re-examination 
of a cascade decay (Bayneres Report, p. 99). In this event, the neutral 
V-particle is identified as a A°, and the plane of the neutral V-particle contains 
the decay point of. the charged V-particle. Assuming this to represent a. 
decay of the form Y-+A° +7--+-Q, the Q-value is calculated to be ~65 Mev. 
In addition, a re-examination (Butler, private communication) of the first 
cascade decay (Armenteros et al. 1952), on the basis of the decay mentioned. 


above, shows a @-value ~15-50 Mev, possibly consistent with the event of 


Cowan. Event Y(B)1 does not therefore appear to be of this cascade type. 


Observations on Charged Hyperons in Photographic Emulsions 423 


case, the Q-values would be reduced by ~2 Mev as compared with the 
case for a neutron secondary. The track length of the primary particles 
is in no case sufficiently long to distinguish between the different hyperon 
masses corresponding to these two possible identities of the neutral 
particle. 

Event Y(B)2 does not fit into either of the decay schemes discussed. 
above. If the possibility of a three-body decay be excluded, then this 
event represents either a decay involving a very high energy release, or a 
nuclear capture with subsequent disintegration. The latter would seem 
the most probable explanation, without necessitating the introduction 
of anew type of hyperon. Since the events described above were observed 
at Bristol, Johnston and O’Ceallaigh (1954) have also observed an event 
which seems to represent the nuclear capture of a negative hyperon. 
This event is described in detail in the following paper. 
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SUMMARY 


An event which is interpreted as a nuclear interaction of a charged 
hyperon has been observed in a photographic emulsion. The appearance 
of the interaction is that of a typical co-star, but is caused by a particle of 
protonic mass. Estimates of mass by two independent techniques are 
in agreement with each other, and are consistent with those of charged 
hyperons reported by other workers. 


§ 1. IyTRODUCTION 


DurRInG examination of plates from a stack of 600, stripped emulsions 
exposed during the Sardinian Expedition of Summer 1953, an event has 
been found which presents the following unusual features. 

A particle of g*~7 is emitted from a star of 30 branches of which the 
majority are ‘ black’ or ‘grey’. Because the star is very close to the 
surface of the emulsion, it has been considered unprofitable, for the 
moment, to undertake a more precise classification of this interaction. 
There are also other reasons which are connected with the processing of the 
plates. These will be discussed in the appendix. It is clear, however, 
that a heavy nucleus is involved, and, further, the energy per track is of 
order 100 Mev. 

Thus, the creation of hyperons (Fretter e¢ al. 1953) appears energetically 
possible. The particle under consideration which resembles a proton, 
traverses 3 emulsions, and stops after a range of 7350. That the 
particle was moving with negligible speed at the end of its track, was 
established by observing the distribution of gaps ; this distribution was 
compared with that of protons brought to rest in the same region of the 
emulsion. A photomicrograph of the event is shown in Plate 14. From 
the point of arrest two particles are emitted. One of these is of very short 
range ~4,y. It appears to lie in a plane parallel to the surface of the 
emulsion. No charged particle can be observed in association with its 
end. The other particle is emitted at a very steep angle. It traverses 
three complete emulsions, but cannot be traced into the fourth. It is 
possible therefore that it came to rest very close to the glass backing plate, 


+ Communicated by the Authors. 
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or in the paper tissue separating the emulsion strips. Since the end could 
not be scrutinized, we cannot exclude the possibility of subsequent 
interaction or decay. Its appearance, however, is that of a singly- 
charged particle more massive than a m-meson, and it seems reasonable to 
suppose that it is a proton. The residual range cannot be measured with 
high accuracy, but it is approximately 2000 u (proton energy ~20 Mev). 

Because of the following observations, it is clear that the event cannot 
be interpreted as a case of elastic scattering of a proton of residual range 
2000 4, with recoil of range 4. The tracks are not coplanar, and momen- 
tum cannot be conserved unless we invoke the emission of at least one 
neutral particle. Moreover, as has been remarked, the gap-density of the 
incoming track at the point of supposed scattering differs significantly 
from that of a proton of the postulated residual range. It is concluded 
therefore, that we are dealing with a nuclear interaction similar to those 
commonly produced by 7-mesons arrested in the emulsion. 


§ 2. Mass or InTERAcTING PARTICLE 


Mass estimates of the interacting particle have been made using two 
independent methods, and these agree within the expected statistical 
uncertainty. 

(1) Multiple-Scattering/ Residual Range 

The method of constant sagitta described by the Brussels and Bombay 
Groups (1953) was used with the following modifications. An estimate of 
noise was made by the method of cross products (Menon et al. 1951, 
Moliére 1951). Large unrepresentative signals were treated by the 
substitution method of Goldsack (1953). The range increment tables 
P=1 and z=1-6 kindly supplied by Dr. Peters were used as a basis for 
these measurements. Evidence for the presence of distortion was found 
from the distribution of the signals and the values of the cross products. 
The necessary correction was small, and was applied by the Cosyns— 
Vanderhaege method using measurements on neighbouring steeply 
dipping tracks. A small correction was also made for the influence of dip. 
The corrected mass-estimates were : 


ere 700 
P=10 22801) jo, Me 
ra G00n 
m=16 21001 jog Me 


(2) Mean Gap-Length/ Residual Range 
The conditions of processing the plates which were cleared after fixing, 
were such that the plateau blob-density was low—14 blobs/l00 pu. In 
consequence, the blob-density/range characteristics were unsatisfactory. 
For protons of residual range Rk, 2mm<Rh<15 mm, the blob- 
density/range curve passes through a flat maximum so that it appears to 
be independent of range. The theory of this effect has been investigated 
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by one of us (O’Ceallaigh 1953 a) and will be discussed in a later com- 
munication. In order to estimate ionization in plates such as the present, 
observations must be made of the variation with range of parameters such 
as the photoelectric density or the mean gap-length G (O’Ceallaigh 
1953 a, Menon and O’Ceallaigh 1954). As hitherto described, neither of 
these methods is convenient or rapid, and the method G/R was therefore 
modified as described in the appendix. 

Estimates of G were made in the three plates through which the 
interacting particle passed, and have been plotted on log—log paper as a 
function of residual range (see figure). 


lOO 


50 


MEAN GAP LENGTH 


iE oe O5 | 2 > (Ores 10 Nhe, 8: 
RANGE (MMS) 


Estimate of mass of interacting particle by G/R. 


The interacting particle was emitted from a star in plate 3, passed through 4 and 
came to rest in plate 5. The experimental points are numbered accord- 
ingly, except for those of the interacting particle measured in plate 3 
which are labelled Y with no subscript. The plate numbers and 
distinguishing marks for the calibration particles are to be found on the 
figure. Gin micrometer drum divisions, 11-4=1 p. 


These values of G have been compared with those found for identified 
m-mesons and stopping protons, of which the tracks passed through the 
suitable regions in the plates. There was found evidence for the presence 
of protons showing anomalously low ionization which may have been 
those of faded or ‘ wet’ tracks. For this reason, the calibration protons, 
with one exception, were chosen from amongst those which could be 
traced back to a nuclear interaction in the emulsion. The values obtained 
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from individual tracks are distinguished by symbols listed on the diagram. 
That the selected particles were very unlikely to have been deuterons or 
tritons was established (a) by observation of the associated 6-rays, (b) 
by measurement of blob-density in regions beyond the extremum, and 
finally (c) by mass estimate by range-scattering using Peters’ method. 

The proton points have been fitted visually by a straight line. This 
line will be seen to be a satisfactory fit. Except for points of small residual 
range where the statistical weight is low, the points are based on segments 
of track ~1000 uw, each containing some 500 visible gaps. Since we have 
not yet investigated the problem of the sampling errors of our estimates of 
G, it has seemed pointless to attempt to fit by an analytical method. 

The z-line has been drawn parallel to the proton line in the position 
calculated from the known mass-ratio, and is in good agreement with the 
experimental points. The open circles marked Y are the experimental 
points for the interacting particle. Of the latter, one point corresponding 
to the segment nearest the star indicated a mass value less than that of the 
proton. That this was probably due to its nearness to the surface of the 
plate, is shown by the pair of points in full circle obtained by dividing 
this segment in two. The point of high @ lies in a region <30 » from the 
surface, and is rejected, the other agrees with the remaining observations. 
A hatched line, corresponding to mass 2300 m, has been drawn parallel to 
the proton line, and appears to fit the points in a reasonable manner. 

The position of the first 6 of these may be compared with those of a 
proton stopping in the same plate (upright crosses). Although the points 
in this region are of inferior statistical weight, the agreement between 
them demonstrates that we are dealing with a particle which was at rest, 
or moving with negligible velocity, at the point of interaction. 


§ 3. DiscussION 


The satisfactory agreement between the mass-estimate by two indepen- 
dent methods, makes it reasonable to suppose that the event observed 
is due to the absorption of a negative hyperon in accordance with some 
mechanism such as 


Ver ee eee er 1) 


Taking the mass of Y~ as 2300 m,, the expected value of Q is about 230 Mev. 
The appearance of the interaction is quite consistent with such value of Q. 
Tt is also consistent with that expected on certain hypotheses (Friedlander 
et al. 1953), when a negative K-particle is absorbed in a nucleus. Although 
the interconnecting particle is not a7-meson, it cannot perhaps be excluded 
that it is a negative K-meson. The agreement between the independent 
estimates of mass is such, however, that it does not seem possible to 
accept this identification without assuming unreasonable values for the 
statistical uncertainties of the measurements. 

The present event may be compared with that reported in this issue by 
Friedlander and called Y(B)2 and with one described by Bonetti et al. 
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(1953). The latter event J, has been interpreted as the decay of a 
hyperon in accordance with the scheme 

Bbc ab amy =) tones tgs lise eae ale Ca 
The proton from J, has a range of 1670418. Thus it, and that des- 
cribed by Friedlander, cannot both be cases of the process (2). It is 
possible, however, to interpret all three events as an interaction of Y~ in 
accordance with (1), although as pointed out by Bonetti e¢ al., owing to the 
short range (900) of the primary particle of J,,, they could not exclude 
the possibility that their event represents the capture of a negative 
K -particle. 
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APPENDIX 


Measurement of Mean Gap-Length (G) Improved Technique 


The rapidity of the technique for estimating G described in earlier 
communications has been much improved by the modifications now to be 
described. In essence, this procedure has been suggested independently 
by Fowler (1953). An eyepiece is furnished with a set of four spider webs 
placed at suitable distances apart. A convenient choice consists of intervals 
of 0-5, 1:0, 2-0 and 4-0 all measured from the first web. Every visible 
gap is observed and may be classified immediately according to which 
intervals it exceeds. Thus, for each segment of track, we may plot on 
log-linear paper, five points of an integral gap-length distribution, and 
after a simple correction for dip, a straight line may be fitted visually to the 
corrected points. Alternatively, it may be shown that if N, is the number 
of gaps exceeding «,, and NV,, the number exceeding ¢,, then the mean gap- 
length G may be found from the relationship (O’Ceallaigh 1953 b) 

Eo 


Wy exe —(e,—«,)/G. ea ERA A eS 


ey) 


The validity of this relationship is unaffected by changes of processed 
grain-diameter, and the elaborate corrections for depth necessary in 
photoelectric techniques are thus avoided. It is necessary, however, to 
exclude observations made very close to the bounding surfaces of the 
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emulsion, as conditions of development are commonly anomalous in these 
regions. 

In practice, it will be found convenient to estimate G from the slope of 
the straight line fitted as described to the experimental points. The 
validity of the method has been tested by comparison with direct measure- 
ment of G using a filar micrometer. The agreement was satisfactory. The 
principal advantage of the method is its rapidity which is an improvement 
by a factor of about 5 on the direct measurement of G by the filar micro- 
meter. Disadvantages are that some information is lost, and that unlike 
those in a direct estimate of G, the sampling distributions remain to be 
investigated. 

It is intended to examine the application of the technique to the classi- 
fication of the branches of stars in circumstances, such as the present, 
where direct blob/gap density measurements do not yield useful 
information. 
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ABSTRACT 


The angular distributions of the high energy gamma-rays (a mixture 
of 3-7 and 3-9Mev gamma-rays) and of the low energy gamma-rays 
(0-2 Mev) from the 1°B(«, py)!#8C reaction, and the angular correlations 
of the high energy gamma-rays with the emitted protons have been 
measured at five alpha-particle resonances at 1-13, 1-51, 1-64, 1-68 and 
1:83 Mev. These experiments confirm the assignment of spins and 
parities 4-, 3- and 4+ to excited levels of 14N at 12-42, 12:69 and 
12-78 Mev, and spins and parities 3/2~ and 5/2+ to excited levels of ¥C 
at 3-68 and 3-89 Mev. 


$1. INTRODUCTION 


We have made a detailed study of the 1°B(«, py)#C reaction with the 
electrostatic generator of the Cavendish Laboratory. A previous article 
(Shire et al. 1953), hereafter referred to as I, gives an account of some 
of these investigations. The present paper describes experiments on 
the angular distribution of the high energy gamma-rays (a mixture of 
3-7 and 3-9 Mev gamma-rays) and the low energy gamma-rays (0-2 Mev) 
from this reaction, and experiments on the angular correlation of the 
short-range protons with the high energy gamma-rays. We have measured 
the distributions and correlations for five alpha-particle resonances at 
1-13, 1-51, 1-64 and 1-83 Mev. From the results we were able to confirm 
the spins and parities assigned in paper I to levels of C at 3-68 and 


3°89 Mey, and to the excited levels of 14N corresponding to the first three 
alpha-particle resonances. 


§2, EXPERIMENTAL METHOD AND APPARATUS 
We made all measurements with a beam of singly charged helium ions 
produced by the electrostatic generator. We used thin targets of 
separated '°B, about 20 ug em-*, on thin silver or aluminium backing. 


The gamma-rays were detected by means of a sodium-iodide crystal 
and a photomultiplier. 


a Rt a ee 


* Communicated by E. S. Shire. 


Angular Correlations in !°B(«, py)3C Reaction 431 


2.1. Proton-Gamma Angular Correlations and Angular Distributions 
of High Energy Gamma-Rays 

For the proton-gamma angular correlations we used two arrangements : 
in the first the gamma-ray detector could be rotated in the vertical plane 
perpendicular to the incident alpha-particle beam ; in the second it was 
rotated in the horizontal plane containing the incident beam. In both 
arrangements we detected the protons with a proportional counter fixed 
in the horizontal plane at right angles to the incident beam. In most of 
our experiments we turned the gamma-ray detector through an angular 
range from 90° to 180° relative to the proportional counter, but a few 
readings were taken at forward angles. The solid angles of acceptance 
of the proportional counter and of the sodium-iodide crystal were usually 
0-12 and 0-15 steradians respectively, but we sometimes used smaller 
solid angles to measure the gamma-ray distributions. 

We used a conventional circuit with a resolving time of 0-2 psec to 
detect coincidences between the protons and the gamma-rays. In,order 
to record true and accidental coincidences simultaneously, we delayed 
the pulses from one of the two channels by 2-5 usec before passing them 
into a second coincidence unit together with the undelayed pulses from 
the other channel (Littauer 1950). The ratio of the resolving times of 
the two coincidence units was determined before and after each 
experiment ; it remained constant over a period of several months. 
We monitored the reaction by means of the proton counts. 

In our experiments we were unable to separate the protons and gamma- 
rays from the following three reactions (see paper I), which we shall refer 
to as (A), (B) and (C): 

(A) 1°B-+ 4He >14N*+8(C* (3-9 Mev) +Pps 
v 
13C (ground state)+- (3-9 Mev). 
(B) 2°B+4He >14N*+8C* (3-9 Mev) +p, 
\ 


180 * (3-7 Mev) ++/(0-2 Mev) 
1 
13C (ground state)-+ (3-7 Mev). 


(C) 1°B-++4He >14N*>130% (3-7 Mev) +p» 


13C (ground state)-+ (3-7 Mev). 
A determination of the branching ratio of the gamma-rays showed that 
the total yields of reactions (A) and (B) were in the ratio 7: 3. The yield 
from reaction (C) was very weak except at the 1-51 Mev resonance, where 
this reaction contributed 16% of the total proton counts. Deuterons 
from the 1°B(«, d)!2C reaction were also detected by the proportional 
counter, but did not produce any true coincidences. 


2.2. Angular Distributions of Low Energy Gamma-Rays 


We studied the low energy gamma-rays by means of an electronic 
pulse analyser of the type described by Hutchinson and Scarrott (1951). 
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We found that the 0-2 mev gamma-rays due to reaction B were super- 
imposed on a large background of low energy pulses. Some of these were 
produced by the generator and the remainder by high energy radiation 
from the target that had been scattered and degraded by the sodium- 
iodide crystal. We reduced the effect of the former as much as possible 
by surrounding the crystal with a thick lead screen, in which we cut a 
hole, 1 in. in diameter, opposite the target. Figure | shows the spectrum 
of gamma-rays of energy less than 700 kev after screening the crystal. 
The two peaks at 200 and at 500 kev were produced by gamma-rays. 
from reaction (B) and by annihilation radiation respectively. 


Fig. 1 
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Spectrum of gamma-rays of energy less than 700 kev at the 1-51 Mev resonance. 
The gamma-rays were detected at 90° to the alpha-particle beam. The 
lower curve shows the background radiation from the generator after 
screening. 

We made measurements with and without a piece of lead foil, 0-115 cm 
thick, placed in front of the aperture. We found that the reduction of 
the low energy background was proportional to the reduction of the 
high energy radiation measured simultaneously. From a calculation 
we estimated that the lead foil should have absorbed 76%, of the incident 
0-2 Mev radiation. From these measurements and from similar measure- 
ments with only the background due to the generator present, we were 
able to determine the amount of 0-2 Mev radiation relative to the amount 
of high energy radiation from the target. 


§3. EXPERIMENTAL RESULTS 
Figures 2, 3 and 4 show the experimental angular correlations and 
gamma-ray distributions after correction for the finite geometry of the 
detectors. We were able to compare the results for the 1-13, 1-51 and 
1-64 Mev resonances with theoretical predictions (see §4); the inter- 
ference produced by overlapping excited states of 14N (see paper I) made 
similar comparisons impossible for the resonances at 1-68 and 1-83 Mev. 


Relative number of true coincidences 
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The angular distributions of the 0-2 Mev gamma-rays are shown in fig. 4. 
The standard deviations of the experimental points are rather large 
for the 1-13 Mev resonance owing to the low reaction yield, and for the 
1-83 Mev resonance owing to the increased background from the generator. 


Fig. 2 
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Proton-gamma angular correlations in plane of incident alpha-particle beam 
(column), in plane perpendicular to incident beam (column 6) and angular 
distributions of high energy gamma-rays (column c). In columns @ and c 
the angles are measured with respect to the incident beam, in column 6 
with respect to the direction of the observed protons in laboratory space. 
————  ----- Theoretical curves. The relative phases of the 
spin-channel amplitudes assumed in reactions (A) and (B) for the 1-64 Mev 
resonance are as shown; for the 1-51 Mev resonance the relative phases 


assumed are: 


Curve Reactions (A) and (B) —- Reaction (C) 
1 0 0 
2 0 7 
3 7 0 
4 7 7 


In b, all the theoretical curves are isotropic within experimental error. 
Magnetic dipole transition has been assumed for the 3-7 Mev gamma-rays 


in all theoretical curves. 


Relative number of true coincidences 
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Angle 
Proton-gamma angular correlations in plane of incident alpha-particle beam 
(column a), in plane perpendicular to incident beam (column b) and 
angular distributions of high energy gamma-rays (columnc). The angles 
are defined as in fig. 2. The smooth curves in column c have been fitted 
to the experimental results by the method of least squares. 


Fig. 4 
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Angular distributions of low energy (0-2 Mev) gamma-rays. 
Theoretical curves, — — —— — — experimental curves obtained by 


method of least squares. The relative phases of the spin-channel ampli- 
tudes assumed for the resonances at 1:51 and 1-64 Mev are as shown. 
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§4. Comparison or EXPERIMENTAL RESULTS AND THEORETICAL 
PREDICTIONS 
4.1. Theory 
Biedenharn et al. (1951) have derived a formula for three-stage angular 
correlations, but this does not cover the case where the emitted particle 
has non-zero spin. We have made use of a method, employed by Spiers 
(1950) for two-stage correlations, in which the angular momentum wave 
function for the final state of the whole system is obtained by successive 
transformations from the initial state. In the final state all the emitted 
particles are at a very large distance from the residual nucleus, so that 
the radial part of the wave function may be omitted from the correlation 
formula. Taking the direction of the incident alpha-particle beam as the 
axis of reference and assuming that only one level of the compound 
nucleus, '*N, contributes to the reaction, the theoretical angular correlation 
between protons emitted in the direction (9,, ¢;) and gamma-rays emitted 
in the direction (@,, dg) is given by 
W(91, $1392, ba)= XX — | Xi"(2) +1)? (SMI, 0 | JM) 
My, M,",M. lo 


EAT, Soy Sas lop Ly) FIM | SM aly) 


can 
x Y"(61, $1) SM, \s, “MS, eaes 
eee riekis Nes its )A (ja) ®;.""(O2, $2) |?. et) 


So is the resultant spin of alpha-particle and target nucleus, !°B, and is 
thus always 3. J is the total angular momentum of the compound 
nucleus, 14N. The spins S,’ and S,” of the excited level of °C and of the 
proton are combined vectorially to produce the channel spin S,. The 
ground level of #C has spin S,=3. The symbols M denote the 
components in the direction of the incident beam of the spins with the 
corresponding suffixes. 


(SoM ol) |JM), (JM |S Mylym,) and (S,M,|S1'M,'S,'M,") 


are Clebsch-—Gordon transformation coefficients, and Y,(0,,¢,) is a 
spherical harmonic for the emitted proton. D(S,’, My’, S,, M,) is a 
transformation coefficient for the multipole field of the gamma-ray and 
®,2""(0, 2) its angular dependence (Condon and Shortley 1951). 
The factor A(J, So, 81, Uo, 11) is the product of the following terms : 
A(J, So; Si, lo, L)=K(, Ly )f (lo, Ey)NU, So, Si, lo, Ly) f( ly, E,) 
(2) 
K(J, E,) is a function of the alpha-particle energy B, wena: expressed 
in the form of a complex resonance denominator. (lp E,) and f(l,, £3) 
are complex factors describing the probability ee Bele 
with orbital momenta, J, and 1,, of the incident and emitted particles. 
The intrinsic nuclear factor N(J, So, 84, 9, ,) depends upon the spins 
and orbital momenta in a wholly unknown way. The dispersion theory 
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of Wigner and Eisenbud (1947) suggests that it is real, but it may be 
either positive or negative. Since the proton has spin 3, the channel 
spin S, may in general take two values, which in turn give rise to two 
values of N(J, So. Sz, lg, 1). In paper I the ratio of the squares of the 
latter has been defined as the spin-channel factor. It may be seen from 
eqn. (1) that the two spin channels interfere coherently, i.e. either by 
addition or by subtraction of the spin-channel amplitudes. The nuclear 
matrix element A(j,) associated with the multipole j, (see Lloyd 1951) is 
independent of the magnetic substates. 

From eqn. (1) we can derive an expression for the angular distribution 
of the gamma-rays by integrating the spherical harmonic for the proton, 
which is not observed, over all space. The resulting angular distribution 
is given by 


WO)= SF Z| Zh’(Mpt+1)¥24S,Mol0 | IM) Z AT, So, Sx, Ip, hy) 


My, My",Ms lym Io S; 

x (JM | 8S, M lym) (SM, | 8y'My'S,"M,") 

x2 D(8y', My’, Sx My)A(5 9) G5" 5) [ee a poe ae 
Je 


In reaction (B) (§ 2.1) two gamma-rays are emitted in cascade. In the 
experiments we observed the angular distribution of the 3-7 Mev gamma- 
rays and their angular correlation with the protons. We have extended 
the theory to a four-stage process, so as to be able to deal with this case. 
The emission of the first (0-2 Mev) gamma-ray, which is not observed, 
leaves the residual nucleus, C, in an excited level with spin S,. After 
the emission of the second (3-7 Mev) gamma-ray with multipolarity j, 
the residual nucleus is in the ground state with spin S,=}. The angular 
correlation between protons emitted in the direction (6,,¢,) and 
gamma-rays emitted in the direction (63, 45) relative to the direction 
of the incident beam is given by 
W(01, $1593, ¢3)= & & | Lt%(21,4+1)/?2(8,Mol,0 | JM) 


Mo, My”, Ms jo,mM2 lo 


x ZA, Sp, Sy lp, 1) Z (IM | S,M1,m,) 


Sy hi M, 
x Y,™(0,, 61) (SM, | S,'M,'S,"M,") 
x D(8,', My’, Sy, M,)A(j) 2 D(Sq, My, Sq, My) 
Is 


LA(JDOP OL GS |E cence een ay 
and the angular distribution of the gamma-rays by 
W(O)= SF EZ | Téh(2y+-1)¥2(SyMyh.0 | JM) 
l 


Mo, My", Ms lh, my jo,mMe Io 
xZ AI, Sq, Sy, Ip, 13) (JM| SM lym.) 
Sy 


x (SM, | 81'S," ," D(Sy, M,', Sy, My) 
x A(jg) 2 D(S., M,, Ss M3)A3( js) P;.(43) [?. . . . . . (5) 
Is 
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4.2. Calculation of Angular Correlations and Distributions 

In order to calculate the angular correlations and distributions we have 
made use of the spins of the excited levels of 14N and 130 and the 
spin-channel factors obtained from the results of the proton distributions 
described in I. We have assumed further that only the lowest values of 
the orbital momenta of the incident and emitted particles consistent 
with the conservation of spin and parity have contributed to the reaction. 
The assignment of multipoles for the gamma-rays is discussed in § 4.3. 
We have combined the angular correlations, and also the angular 
distributions of the high energy gamma-rays in reactions (A), (B) and (C) by 
means of weighting factors derived from the relative yields of the 
different groups of protons and gamma-rays. The weighting factors 
were corrected for the anisotropic emission of the different particles. 
The results of the calculations are shown in table 1. 


Table 1. Theoretical Angular Correlations and Gamma-Ray Distributions 


Relative Angular Correlation in Angular Distribution of 
i esonance Phases of Plane High Energy 
coeee Spin-Channel Plane containing Perpendicular|Gamma-Rays| Low Energy 
: BY Amplitudes Incident Beam to Incident | (Mixture of |Gamma-Rays 
ak in Reactions Beam 3:7 and 3:9 (0-2 Mev) 
(A) & (B)| (C) Mev) 
1-13 — — 1+0-037 cos (20+21°) | 1+0-24 cos*d | 1+0-17 cos? | 1— 0-27 cos?0 
+0:012 cos (486—17°) 
1-51 0 0 1+0-088 cos (26—78°) | 1+0-01 cos?¢ | 1—0-24 cos?@ | 1-++0-27 cos?0 
| +0-088 cos (40—18°) 
0 7 1+0-118 cos (28—47°) | 1—0-02 cos?¢ | 1—0-19 cos?6 
+0-088 cos (48—18°) 
7 0 1—0-133 cos (26+15°) | 1+0-02 cos?¢ | 1—0-35 cos?6| 1+0-52 cos? 
+0-023 cos (48—18°) 
| 7 7 1—0-074 cos (26+28°) | 1—0-01 cos*d | 1—0-31 cos?6 
+0-023 cos (49—18°) 
1:64 0 — 1+0-064 cos (26+20°) | 1—0-26 cos? | 1+ 0-19 cos?@ | 1 —0-02 cos?0 
—0-026 cos (49+25°) | +0-08 costd | —0-27 cost 
7 — 1—0-095 cos (26+17°) | 1—0-17 cos? | 1+ 0-25 cos?6 | 1—0-03 cos*d 
+0-050 cos (48—10°) | +0-05 cos*f | —0-20 cos*é 


Magnetic dipole transition has been assumed for the 3-7 Mev gamma-rays 


If we take (7/2, 0) as the direction of the protons in eqn. (1), the 
spherical harmonic Y,”"(7/2, 0) causes many of the cross terms to cancel, 
so that the resulting expression becomes particularly simple. But if 
the protons are observed at 90° to the incident beam, they will be moving 
in a backward direction in centre-of-mass coordinates, and then the 
correlation curve in the plane containing the incident beam may be 
substantially different from the curve for which the protons leave at 90° 
in centre-of-mass coordinates. If the correlation is given by 

W(0,)=1-La cos? 0,+6 cos* 0, (6) 
when 6,—90°, then the change of 6, to 90°+-« (where « is the centre-of- 
mass correction to 6,) introduces two additional terms csin 20, and 
d sin 40,, as well as corrections to a and b which are small enough to be 
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ignored. On rearranging the terms, the corrected correlation becomes 
b\2 1/8 
W(6,)=14+4b+3a+ (>) +0 cos (20.—y) 
2 1/2 
+ | (5) +e] cos (40,—a) > 8 


where tan v=8d/a ; tan y=2c/(a+b). Ifa and b are approximately equal 
and of opposite sign, even a small value of ¢ may alter the shape of the 
correlation curve considerably. The effect of the proton motion was 
found to be particularly important for reaction (A) at the 1-51 Mev 
resonance. 
4.3. Confirmation of Spins and Parities 

In paper I it was shown from the proton angular distribution experi- 
ments and from the existence of a strong 0-2 Mev gamma-ray that the 
3-68 and 3-89 Mev levels in %C.must have spins and parities 3/2~ and 5/2+ 
respectively. The observed angular distributions of the 0-2 Mev 
gamma-rays indicate that this gamma-ray is due to an electric dipole 
transition, consistent with the above assignments of spins and parities. 
The angular correlations ana distributions of the high energy gamma-rays. 
showed terms up to the fourth power of cos @, agreeing with the assign- 
ment of 5/2+ for the level at 3-89 Mev; gamma-rays emitted from a 
level with lower spin cannot give rise to terms higher than cos? 0@. 


Table 2. Comparison of Experimental and Theoretical Angular 
Distributions of High Energy Gamma-Rays at 1-51 Mev Resonance 


Experimental distribution 1—(0-151 +0-007) cos? 6 
Theoretical distribution assuming magnetic 
dipole transition for 3-7 Mev gamma-rays 1—0-19 cos?6 
Theoretical distribution assuming electric 
quadrupole transition for 3-7 Mev gamma- 
rays 1—0-29 cos? 6 


The experimental distribution has been obtained by the method of least 
squares. The theoretical distributions are in closest agreement with the 
experimental distribution for which the relative phase ofthe spin-channel 
amplitudes has been assumed to be 0 in reactions (A) and (B), and 7 in 
reaction (C). 


The transition from the level at 3-68 Mev to the ground state might be 
by mixed magnetic dipole and electric quadrupole radiation, but the 
relative transition probabilities based on the independent particle model 
(Blatt and Weisskopf 1952) favour the magnetic dipole transition by 
400: 1. In the theoretical calculations we have assumed a pure magnetic 
dipole transition for the 3-7 Mev gamma-rays, but at the 1-51 Mev 
resonance we have also considered a pure electric quadrupole transition. 
The theoretical expression based on the magnetic dipole transition 
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produced closer agreement with the experimental curve for the gamma-ray 
angular distributions (see table 2), but the influence of the 3-7 Mev 
gamma-rays on the distribution was much reduced by the preponderance 
of the 3-9 Mev gamma-rays. 

Our confirmation of the assignments made in paper I for the spins and 
parities of the excited levels of 14N is only indirect and depends on the 
reasonably good agreement between the theoretical curves derived by 
assuming these assignments and the experimental results here presented. 


4.4. Assignment of Phases 

In §4.1 we have assumed that the amplitudes of the competing spin- 
channels are either in phase or 180° out of phase. It is in principle 
possible to determine the relative phases by comparing the theoretical 
and experimental curves. At the 1-51 Mev resonance the spin-channel 
amplitudes for the predominant reactions (A) and (B) are in phase. This 
is clearly shown by the angular distribution of the high energy gamma-rays 
and by the correlation in the plane containing the incident beam, fig. 2 ; 
it is confirmed by the angular distribution of the 0-2 Mev gamma-rays, 
fig. 4. The distribution of the high energy gamma-rays also suggests 
that in reaction (C) the spin-channel amplitudes are 180° out of phase ; 
but we cannot make a unique assignment, as this reaction produces 
only 16% of the total number of counts. At the 1-64 Mev resonance the 
angular distribution of the high energy gamma-rays indicates that the 
spin-channel amplitudes in reactions (A) and (B) are 180° out of phase. 


$5, CONCLUSION 


Our results confirm the assignments made in paper I for the spins 
and parities of excited levels of ‘4N and #C. The amplitudes and phases 
associated with the channel spins resulting from the spin orientation of 
the emitted proton have been calculated by K. V. Roberts (private 
communication) from the nuclear shell model with Russell—Saunders 
coupling (see Christy 1953), but agreement with our experimental values 
was poor. Recently, Lane (1953) has shown that the low excited states 
of 3C are best explained by means of a nuclear shell model with 
intermediate coupling between the Russell-Saunders and j-j extremes 


(see also Inglis 1953). 
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Photomicrographs (a), (b), and (c) show the V-events (1), (2), and (3) of table 1. 
The values QYg. are the Q-values given in col. 10 of Table 1; the values 
(J,« are those of col. 6 of table 2. 
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Event Y(B)1 : Kijected from a star of type 26+14p, the hyperon decays in 
flight after travelling a distance of 4:72 mm. 
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Event Y(B)2: The hyperon comes from a star of type 31+3n, and reaches 
the end of its range after a distance of 4:64 mm. From the end of the 
hyperon track, there starts another track which can be identified as that 
of a proton, of energy 100415 Mev, 
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